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ABSTRACT 

We introduce the Swift Gamma-Ray Burst Host Galaxy Legacy Survey (“SHOALS”), a multi¬ 
observatory high-redshift galaxy survey targeting the largest unbiased sample of long-duration gamma- 
ray burst hosts yet assembled (119 in total). We describe the motivations of the survey and the 
development of our selection criteria, including an assessment of the impact of various observability 
metrics on the success rate of afterglow-based redshift measurement. We briefly outline our host- 
galaxy observational program, consisting of deep Spitzer/IRAC imaging of every field supplemented by 
similarly-deep, multi-color optical/NIR photometry, plus spectroscopy of events without pre-existing 
redshifts. Our optimized selection cuts combined with host-galaxy follow-up have so far enabled 
redshift measurements for 110 targets (92%) and placed upper limits on all but one of the remainder. 
About 20% of GRBs in the sample are heavily dust-obscured, and at most 2% originate from z > 5.5. 
Using this sample we estimate the redshift-dependent GRB rate density, showing it to peak at z ~ 2.5 
and fall by at least an order of magnitude towards low (z = 0) redshift, while declining more gradually 
towards high {z ~ 7) redshift. This behavior is consistent with a progenitor whose formation efficiency 
varies modestly over cosmic history. Our survey will permit the most detailed examination to date of 
the connection between the GRB host population and general star-forming galaxies, directly measure 
evolution in the host population over cosmic time and discern its causes, and provide new constraints 
on the fraction of cosmic star-formation occurring in undetectable galaxies at all redshifts. 

Subject headings: gamma-ray burst: general — galaxies: star-formation — galaxies: evolution — 
galaxies: high-redshift — surveys 
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Long-duration gamma-ray bursts (GRBs) are rela¬ 
tivistic, jetted explosions of very massive stars at the 
end of their lives.The peak luminosities of GRBs 
and their afterglows exceed those of the most luminous 
supernovae , galaxies and qua sars by orders of magni¬ 
tude (e.g., iBloom et al.l l2009l l throughout the electro¬ 
magnetic spectrum. As a result GRBs are routinely 
detected from great distances: the median redshift of 
Swift G RBs is z ^ 2 (|,Iakobsson et al.ir200(il : iFvnbo et al.l 
l2009bll and G RBs have been de t ected up to a redshift 
of z = 8 — 9 (Tfanvir et al.l 120091 : ISalvaterra et al.l [20091 : 
iGucchiara et al.ll2011h . 

GRBs are therefore fundamentally cosmological 
objects, and their study is intimately coupled with that 
of high-redshift galaxies and cosmic evolution. GRB 
afterglows make excellent probes of individual galaxy 
sightlines, a technique that has been extensively ex¬ 
ploited to characterize the interstellar medium in distant 
star-forming galaxies—chemical abundances, kinemat¬ 
ics, and dust properties—to a le vel of detail unmatched 


20091: lEliasdottir et al. 

120091: iD’Elial 20111 Zafar et al. 

2012 iDe Gia et al.ll201 

IllThdne et al.ll201.‘lHSDarre et al. 


UJ). They 

may also be useful for studies of galaxy and cosmic 
evolution in a broader sense. Because GRBs originate 
from short-lived massive stars, they stochastically 

They are physically distinct from short (Tgo < 2 sec) gamma- 
ray bursts, w hich are belie ved to be associated with merging com¬ 
pact objects IIBergerfl2014l') and are not discussed in this work. 
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sample the sites of cosmic star formation in proportion 
to their relative contributions to the cosmic total, 
providing (in principle) a means of estimating the 
importance of different epochs or different galaxy pop¬ 
ulations to th e Universe’s s t ellar m ass assembly (e.g., 
Totanil Il997t IWiierset^ Il998i iBlain fc Nataraia^ 


200Ct iPorciani fc Madaul l2001li . This includes popu¬ 


lations difficult to s tudy by othe r means: very low- 
lumin o sitv galaxies (iBerger et al.l l2002t iSchulze et alJ 
1201 5t rdreiner et alJ ^1^. very dusty ga laxies 


( Diorgovski et alJ 12001 : iRamirez-Ruiz et al.l 120021: 


20061 

Michalowski et al.l 

20081: iPerlev et al.ll2009cl 1201.31 

2015a 

Rossi et al. 120121: 

Greiner et al. 

201lUHunt et al.l 

2011 

Schadv et al.l 2014 

: iKohn et al. 

201511. verv high- 


l2012t iTrenti et alJ 201211 . and possibly even Populat ion 


111 stars ( Bromm &: Loebir2006i ICampisi et ^l2011a|l . 


These broader questions of using GRBs as probes of 
galaxy evolution and cosmology have provided major 
drivers for GRB research over the past two decades and in 
part icular since 2004, wh en the launch of the Swift satel¬ 
lite (jGehrels et al.l[2003 l revolutionized the field with its 
ability to detect and instantly localize GRBs in large 
numbers (^^100 per year). The vast (and still-expanding) 
Swift legacy data-set, including follow-up provided by 
the world-wide ground-based community, has made pos¬ 
sible the study of GRBs, their environments, and their 
host galaxies in sufficiently large numbers to investigate 
these topics statistically. 

The cosmic star-formation rate is not necessarily 
the only factor expected to affect the GRB rate, 
however. The initial chemical composition of a 
star influences its fate, and many models specifi¬ 
cally predict that the GRB rate in metal-poor envi¬ 
ronments should be diffe rent—typically, higher—than 
in metal-rich ones (e.g.. iMacFadven fc Woosl^ 119991: 
Fryer fc Hegeil 120051: iHirschi et al.l 120051 : lYoon fc Langed 


20051 : IWooslev fc Heged l2006ll . Alternatively, varia¬ 


tions in the IMF o r close-binary fraction could af - 
fect the GRB rate (|Davd 120081 IWang fc Pail I2011I1 . 
as might dynamical interactions in d ense clusters 
(jvan den Heuvel fc Portcgies Zwai^ 1201,811 . If present, 
these effects would have to be taken into account when 
attempting to employ GRBs to investigate broader cos¬ 
mological issues. 

A non-constant GRB production efficiency (defined 
as the GRB rate relative to SFR) is supported by 
several studies of the GRB host population at z < 
1 showing known GRB hosts to be smaller, bluer, 
less massive, and less chemically enriched on average 
than typical star-forming galaxies or core-collapse super¬ 
nova hosts (iFruchter et a.l1 l2006t ILe Floc’h et al.l 120061: 
Stanek et al] 120061: iModiaz et al.l 12 0081: ISvensson e t al l 


20101: iLevesoue et a,l.l 120101: iGraham fc Fruchted 1201.81 

Boissier et al.l 120131 : iKellv et al.l I20l4 IVergani et al.l 

201, 5fl . Furthermore, a number of studies have concluded 


that the comoving GRB rate density (as inferred from the 
redshift distribution) evolves differently from the star- 
formation rate density, showing an excess of GRBs orig¬ 
inating from z > 3 compared to what would be expected 
if the GRB rate was a fixed fraction of the star-formation 
rate in all types of galaxy and at a l l redshifts (e.g., 
iDaigne et al.l l2005 iLe fc Permed l200'^ iGuetta fc PiranI 


20071: iSalvaterra & Chincarinil 120071: lYiiksel et al.l 120081: 

Kistler et al.ll2008ll. 

re connection be- 
ation rate is not 
n whether the na- 
>t explained by a 
)ther factor which 
nation rate inten- 
al.l 120141. Other 
rate varies sig- 
h redshift at all 

Even so, the actual nature of tl 
tween the GRB rate and star-form 
yet well-determined. It is not certai 
ture of the low-z population is bes 
dependence on metallicity or on anc 
correlates with it, such as star-fom 
sitv (iKocevski & West! 120111: IKellv et 
authors dispute whether the GRB 
nificantlv with host galaxv or wit 

le.g.. IWanderman & PiranI 120101: iCai 

mnisi et al.l 1201 Ibl: 


Elliott et al.ll2011 iMichalowski et al. 

120121 iHunt et al.1 

2014IKohn et al.N201.5ll. 

rs that affect these 
rests are those of 
hat dominate the 
Lve properties very 
rergetics, intrinsic 
the cosmological 
[shifts (iPian et al.l 

There are several complicating facto 
studies. Many of the best-studied 
the so-called low-luminosity GRBs 1 
rate locally (z < 0.25)—but which ha 
different (orders of magnitude in ei 
rate, and degree of beaming) from 
GRBs that are detected at high rec 

20061: IGuetta & Della Valid 120071 iWiersema et al.l 120071 


Virgili et al.l 120091 : iBromberg et al.ll2011ll . Despite their 


much greater observed numbers, the hosts of the dis¬ 
tant, cosmological GRBs are more difficult to study and 
have been explored less consistently and in much less 
detail, and major selection biases underly our ability to 
produce representative samples. The studies of the cos¬ 
mological GRB population that have been conducted to 
date have in all cases been small (< 20 objects), heav¬ 
ily selection-biased, or characterized the hosts only min¬ 
imally (observations in 1-2 photometric bands, often to 
limited depth.) 

In this paper we introduce and summarize the largest 
and most ambitious survey of the GRB host galaxy pop¬ 
ulation conducted to date, designed to move past these 
earlier limitations by constructing a sample which is 
large, unbiased, and thoroughly observed at a variety of 
wavelengths across the electromagnetic spectrum. Our 
effort, which we designate the Swift gamma-ray burst 
HOst gAlaxy Legacy Survey (or “SHOALS”), provides 
the most complete view (in wavelength, redshift, and 
sample construction) yet of the environments in which 
GRBs explode in order to rigorously examine all aspects 
regarding the connection between gamma-ray bursts and 
galaxies across cosmic history. What is the GRB rate 
history? What is the distribution of properties among 
galaxies that host GRBs, how does it compare to the 
properties of galaxies thought to dominate the Universe’s 
star-formation activity, and how does this change with 
redshift? What fraction of the population originates 
from obscured regions or from dusty galaxies? What 
fraction originates from galaxies too faint to be detected 
at all? What can we learn by combining properties 
measured by observing GRB afterglows in absorption 
(dust columns, kinematics, and gas/metal columns) with 
those learned by studying their hosts in emission (mass, 
star-formation rate, morphology, and nebular abundance 
measurements)? Our survey seeks to answer these and 
similar questions. 

This paper provides an introduction to the survey, in¬ 
cluding an explanation of our selection methodology ((J2]~ 
IS), a summary of our observational campaign ((j4]), pre¬ 
sentation of our new redshift determinations ((J5]), as well 
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as tables detailing the properties of the sample. Using 
these observations, which establish the largest highly- 
complete GRB redshift sample to date, we provide a 
measurement of the inferred cosmic GRB rate density 
(ijH) unbiased by selection effects, and discuss its rela¬ 
tion to the cosmic star-formation rate density (IJT]). We 
summarize these results and the project in general in iJS) 

Subsequent papers will summarize science results asso¬ 
ciated with the host-galaxy campaign, an ongoing effort 
to provide complete SED measurements for all our tar¬ 
gets. In particular, our Spitzer observations and the NIR 
luminosity and stellar mass distribution inferre d from 
will be presented in Paper II (|Perlev et al.ll201^ . which 
is submitted concurrently with this work. Additional 
papers will focus on other science questions, including 
the UV luminosity distribution and its evolution as in¬ 
ferred from the optical photometry, careful determina¬ 
tion of the host-galaxy physical properties (SFR, mass, 
etc.) from SED fitting and an assessment of their im¬ 
pact on the GRB rate, detailed analysis of correlations 
between properties of a host and properties of its after¬ 
glow, the spectroscopic properties of the host galaxies, 
an investigation into candidate hosts of foreground ab¬ 
sorbing systems seen along GRB sightlines, and further 
topics. 

2. SAMPLE CONSIDERATIONS 
2.1. Biases Affecting Redshift Measurement 

Swift detected 803 long-duration (Tgo > 2s) GRBs 
from the start of the mission through to the end of 2014, 
and approximately 320 of these events have measured 
redshifts: a potentially enormous sample to draw from. 
However, from the point of view of characterizing the 
host population or redshift distribution, attention must 
be paid to systematic biases as well as raw number statis¬ 
tics. In particular, the subset of Swift GRBs with mea¬ 
sured redshift is not only incomplete (40% of the popula¬ 
tion) but very likely to be biased: the ease of providing a 
redshift measurement for a GRB is expected to be depen¬ 
dent on properties of its h ost galaxy and on the redshift 
itself (jCoward et al.l[2M^ . 

The vast majority of GRB redshifts in the Swift era 
are provided by absorption spectroscopy of the opti¬ 
cal afterglow, and GRB afterglows show great variety 
in their optical brightnesses, especially at early times 
(|Akerlof fc SwanI l2007t iKann et ^ 1201011 . In particu¬ 
lar, the afterglows of some events are sufficiently faint 
{R > 22 mag, even during the first hour) that they 
cannot be detected by the small-to-moderate aperture 
telescopes typically used to identify the optical after¬ 
glow before triggering large spectrographs (and in some 
cases, even deep optical imaging at large telescopes fails 
to detect an afterglow). These are typical l y referred to 
as “d ark” GRBs (e.g., iGroot et al.lIl99d iFvnbo et al.l 
l2001ll ^°: redshift catalogs based primarily on publicly- 
reported afterglow redshifts necessarily exclude them and 

Alternatively, a GRB can be defined as “dark” based on the 
degre e of optical faintness relative to X-ray waveleng ths (e.g. via 
/3ox; IJakobsson et aO2004l : Ivan der Horst et al.ll200^ — which iso¬ 
lates optically-absorbed bursts specifically, including very luminous 
afterglows that can shine through thick dust screens but excluding 
GRBs whose faintness is intrinsic. Since this work is focused on 
afterglow sample completeness (and not necessarily its underlying 
causes) the former, “absolute” definition is generally more relevant 
here. We employ the term loosely to refer to GRBs with optical 


are therefore intrinsically optically-biased. This bias re¬ 
moves GRBs in dusty environments or at very high red¬ 
shifts, two of the most interesting regimes we might wish 
to use GRBs to explore, and without which our view of 
the redshift distribution and host population are severely 
incomplete. 

Fortunately there are other means available to us to 
recover the GRB redshifts even without an optical de¬ 
tection: GRB afterglows are luminous across the entire 
electromagnetic spectrum, and the afterglow position can 
often be recovered via observations at near-infrared, X- 
ray, or radio wavelengths, permitting the host galaxy to 
be identified and its redshift measured spectroscopically 
in emission (or, less optimally, photometrically). Histor¬ 
ically, subarcsecond follow-up at wavelengths outside the 
optical band was much less routine than optical follow-up 
and the problem was nearly intractable from a statistical 
point of view: dark GRB hosts and r edshifts could be 
recov ered in a few special cases (e.g., iDiorgovski et al.l 
I200II1 but most of them were lost forever without an ac¬ 
curate localization. 

This situation was transformed by Swift. The satel¬ 
lite is equipp ed with an on-board focusing X-Ray Tele¬ 
scope (XRT: iBurrows et al.ll200dl and the capability of 
automatically slewing to a field immediately after de- 
tecting an event with its Burst Alert Telescope (BAT; 
iBarthelmv et al.l l2005li , meaning that early time, high 
spatial-resolution X-ray coverage is available for every 
event the BAT detects. XRT observations of Swift GRBs 
almost always detect the afterglow: 98% of long GRBs 
are detected by XRT (iBurrows et al.l[200^ . and if events 
for which the XRT follow-up was delayed are excluded 
this figure rises to 100%. So even without an optical 
counterpart it is almost always possible—with sufficient 
effort—to recover (or place a deep upper limit on) the 
host galaxy using the ~ 1.5" — 2.0" X-ray positions rou¬ 
tinely delivered by the XRT (jButlerl [20071 iGoad et al.l 
120071: lEvans et ahl 120091) . These rapidly-available posi¬ 
tions also make it much easier to identify faint optical 
counterparts (or to recognize their absence and trigger 
deeper optical follow-up or multi-wavelength observa¬ 
tions) for bursts whose “darkness” is borderline. 

Host-galaxy spectroscopy (and host-galaxy observa¬ 
tions in general) is observationally expensive and can 
only be carried out for a limited number of objects: re¬ 
covery of all missing Swift redshifts is unfeasible. More¬ 
over, absorption redshifts can be missed for a variety of 
reasons other than intrinsic darkness: GRBs occur at 
random times and directions and some cannot be fol¬ 
lowed if, for example, the location is too close in projec¬ 
tion to the Sun or the Galactic plane. Indeed, about half 
of Swift GRBs are probably missed for mundane reasons. 
Several studies have estimated an intrinsic dark fraction 
of about 20% (ICenko et al.l l2009bl : iGreiner et al.l [ 2 M 1 I : 
iMelandri et al.ll2012D . so 80% of Swift GRBs should have 
afterglows bright enough in principle for optical spec¬ 
troscopy, yet in only 30% of cases are redshifts actually 
reported. Expending effort and observational resources 
to recover hosts and redshifts for events missed for these 
reasons would not be informative. 

At the same time, for the reasons discussed above, 

afterglows that are sufficiently faint that afterglow-based redshift 
recovery is unusually challenging or impossible. 
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we cannot simply neglect GRBs without afterglow red- 
shifts without also discarding obscured or high-z events 
from the sample. Two general approaches are available 
to characterize the complete population. The first is to 
simply take the known-z population and attempt to “cor¬ 
rect” it by also observing a limited number of GRBs 
that are definitely dark (have deep, early optical lim- 


in the past (e.g., 

Kriihler et al.l 120111: iRossi et al.l 120121: 

iPerlev et al.l 12013 

: Hunt et al.l 1201411 and among other 


things has clearly shown that GRBs with very faint op¬ 
tical afterglows probe a different (dustier and more mas¬ 
sive) host population from GRBs with bright optical af¬ 
terglows. However the actual contribution of dark GRBs 
to the whole is not precisely determined, in part because 
the matter is much more complicated than an either/or 
distinction between “bright” and “dark” bursts: the dis¬ 
tribution of afterglow brightnesses is a continuum, as is 
(we expect) the likelihood of redshift measurement as a 
function of brightness. “Very” versus “somewhat” dark 
bursts may likewise have different host and redshift dis¬ 
tributions. 

An alternative approach, and the one we adopt in this 
work, is to carefully down-select the Swift sample to 
remove GRBs that occurred under circumstances that 
were not optimal for ground-based follow-up and iso¬ 
late a sub-set for which the afterglow redshift complete¬ 
ness is close to the expected maximum achievable value 
of about 80% (the remaining 20% being dark bursts). 
This basic technique was first exploited to study unbi- 
ased afterglow demographics and redshift distributions 
(iJakobsson et al.ll2006t iFvnbo et al]l2009bt ICenko et al.l 


l2009bHPerlev et al.l2009cl : lGreiner et al.ll2011h . and more 

recently has been successful in addressing the proper¬ 
ties of the GRB host galaxy population as well. Most 
notably, the multi-year VLT-based “TOU GH” project 
(The Optically Unbiased GRB Host Survey: iHiorth et all 
l2012fl used a series of obseryability cuts to isolate a sam¬ 
ple of 69 objects out of the broader Swift sample and 
has achieved a redshift completeness of close to 90% af¬ 
ter intensive optical spectroscopy of the host galaxies of 
those events lacking bright afterglows JJakobsson et al.l 
l20Tl[K7iihler et al.ll2nia ISchulze et al.lf2015ll . “BAT-6” 
(jSalvaterra et al.1l2012tl is an effort with a similar size, 
design, and redshift completeness (58 objects at >90% 
completeness.) 

Though these projects have been informative in char- 


actering the host p o pulati o n in various ways (e.g. , 
Milvang-Jensen et al.l 120121: iMichalowski et al.l 120121 
Schulze et al.ll2015l IVergani et al.ll2015l l and have pro¬ 
vided the firs t nearly-unbiased estim ates of the redshift 
distribution (IJakobsson et al.ll2012f) . from the perspec¬ 
tive of measuring actual redshift evolution of the GRB 
host population, TOUGH and BAT-6 are limited by their 
modest sample sizes. For example in the case of TOUGH, 
the 69 events span a redshift range of 0 < z < 6.3 over 
which an enormous degree of cosmic evolution has oc¬ 
curred: to provide a snapshot of the host population at 
a particular epoch, and compare to other epochs, would 
require splitting the sample into at least 5-6 redshift bins, 
each of which would have only a few host galaxies. This 
is enough to make broad statements, but less than nec¬ 
essary to characterize the parameter distribution at any 
point in history. 


Larger samples are therefore needed. In principle this 
could simply be done by extending the year cut-off affect¬ 
ing earlier samples closer to the present (e.g. TOUGH 
was limited to events from 2005-2007 based on the time 
the survey was conducted; this cutoff could be extended 
to include more recent bursts). However, we can now 
do better: while the detailed selection criteria for earlier 
samples were developed based on well-informed guessing 
about the factors that affect the success rate of after¬ 
glow follow-up, we now have the benefit of hindsight in 
the form of a large Swift GRB population from which to 
establish in practice what parameters and values maxi¬ 
mize the benefit to our science goals. 

As a result, before proposing for and executing our 
survey we devoted significant attention to examining the 
impact of a variety of criteria on the size and complete¬ 
ness level of the resulting sample, deciding on a final set 
of optimized parameters based on these investigations. 
As the associated findings are relevant to understanding 
the design of our survey and may also be useful to inves¬ 
tigators considering even more ambitious projects in the 
future, we describe these in detail in the next section. 

2.2. Additional Considerations 

Since it is essential to produce an unbiased popula¬ 
tion, we consider only observational criteria that are not 
expected to be physically connected to the burst’s en¬ 
vironment in any way. Given the cosmological nature 
of GRBs, factors related to time or sky location are 
all clearly unconnected to the GRB host environment 
and therefore fair to consider. Properties intrinsic to 
the GRB itself (measurements of the prompt emission 
or afterglow) are less obvious. Factors specific to the 
prompt emission should not be biased with respect to en¬ 
vironment, sinc e in the favored in ternal/external shock 
paradigm fe.g.. ISari fc Pira^ll997ll the physical mecha¬ 
nism producing the GRB prompt emission is unrelated 
to the nature of the environment it formed in.^^ On 
the other hand, the GRB afterglow is produced by in¬ 
teraction of t he GRB ejecta with th e circumprogenitor 
environment (|Meszaros fc Ree^ll997f ) and its light may 
be attenuated by gas and dust within its host galaxy, so 
we do not consider any afterglow-related property. 

Because our selection cuts are developed a-posteriori, 
it is important to note that even the application of nom¬ 
inally unbiased cuts can produce some degree of bias in 
a discrete sample, depending on how many features are 
considered and how precisely fine-tuned the cut threshold 
is optimized. If the cuts are drawn from a large potential 
parameter space and optimized directly against the data 
with no separate training set, the choice of cut param¬ 
eters and their values may become driven by stochastic 
effects—artificially driving up the completeness and re¬ 
introducing biases associated with redshift measurement. 
To minimize this risk, we only examine features for which 
a link between the variable and redshift completeness is 
at least potentially expected from fundamental consider¬ 
ations, and only cut on features where significant depen¬ 
dence is observed. We also only consider round-number, 
discrete thresholds in determining the cut value. (An 

This does assume that long GRBs all represent a single class of 
object, rather than a superposition of multiple types with different 
progenitors and prompt-emission properties. 
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alternative would have been to divide the sample into a 
training set and the applied sample, but as the number of 
known-redshift GRBs remains modest—a few hundred— 
this would overly restrict the sample size of both sets.) 

Ghoosing how many cuts to apply, and how strin¬ 
gently to cut, also requires striking a balance between 
redshift completeness and overall statistical size. Even 
when starting with an initial sample of many hundreds 
of Swift events, discrete cuts are necessarily blunt and 
probabilistic tools and many known-redshift events will 
be lost in trying to increase the completeness significantly 
from its very low initial value of 30%. We set as a reason¬ 
able goal obtaining at least 100 events, and at least 65% 
initial (pre-host-follow-up) redshift completeness. The 
target sample size is needed to provide a significant ex¬ 
pansion over earlier host surveys and to ensure that the 
sample is large enough such that, even if subdivided into 
5 or more redshift bins, there are enough targets per bin 
to reasonably statistically constrain the parameter distri¬ 
bution. The completeness goal is informed largely by the 
success of TOUGH in increasing its redshift completeness 
from a similar starting level to its current ^90%, which 
is about the level that is necessary for systematic consid¬ 
erations (preferentially missing redshift measurements of 
faint host galaxies) not to dominate the statistical ones. 
(We had no trouble meeting these goals, with our final 
cuts producing a sample of 119 events at pre-host redshift 
completeness of 68%.) 

Redshift completeness is not the only motive for ap¬ 
plying cuts to the Swift sample—it is also desirable to 
maximize overlap with existing observations of GRBs 
and their hosts. High-quality early afterglow observa¬ 
tions, even if not leading to a successful redshift mea¬ 
surement, can produce a more secure host identification, 
a redshift upper-limit, and useful complementary infor¬ 
mation about the GRB sightline; previous host-galaxy 
observations reduce the observational demands needed 
to complete the survey. 

Because our Spitzer/IRAC campaign f H4.1l) was con¬ 
ducted during Gycle 9 (November 2012 through Septem¬ 
ber 2013), our targets were necessarily restricted to 
bursts occurring before that period, so we only consid¬ 
ered events up to the end of October 2012 for inclusion. 

3. SAMPLE SELECTION AND ANALYSIS OF FACTORS 
INFLUENCING REDSHIFT COMPLETENESS 

With the above considerations in mind, we investigated 
the impact of many different observables to identify what 
combination of unbiased cuts would produce the highest 
final redshift completeness. These are explained in detail 
in the following sub-sections, with the results summa¬ 
rized in Figure [1] Our data were drawn from the Swift 
GRB table^^, supplemented by our own corrections to 
redshifts or classifications where appropriate, using pub¬ 
lished sources and our own observations. 

For the purposes of this analysis we did not treat red- 
shifts that were determined at late times via host-galaxy 
observations as “known”, since the latter subset may be 
affected by conscious (and potentially biased) follow-up 
considerations. We did, however, consider emission-line 
redshifts that were promptly reported in the GCN cir¬ 
culars. (For low-redshift, low-luminosity GRBs the host 

http://swift.gsfc.nasa.gov/archive/grb_table/ 


galaxy may be of comparable or greater brightness than 
the optical afterglow by the time the latter is observed, 
so an emission-based redshift may be reported even if an 
absorption redshift may have been possible in absence of 
host emission lines.) 

We consider only GRBs detected onboard the satel¬ 
lite, excluding ground or slew-survey triggers (which are 
distributed only after many hours’ delay and are rarely 
followed). 

3.1. Time of Explosion 

The amount and quality of observational follow-up di¬ 
rected towards a given GRB afterglow is affected by ter¬ 
restrial timescales in many ways. Observatories are geo¬ 
graphically clustered (most large optical facilities are lo¬ 
cated in Hawaii, the western United States, Chile, or the 
Canary Islands) and the time of day a burst explodes (as 
well as its location) affects which observatories can follow 
it and how quickly. Seasonal weather patterns also influ¬ 
ence the likelihood of whether a given facility is able to 
follow-up a GRB (and, if so, how effectively). On longer 
timescales, other factors may come into play: interest 
in GRB follow-up among different observational groups 
has both waxed and waned over Swift^s lifetime, while 
ground-based facilities important to GRB follow-up have 
been commissioned (e.g. GROND, RATIR, GTC, X- 
shooter) and decommissioned (e.g., PAIRITEL, LOTIS) 
over the same period. 

In Figure we plot the redshift completeness as a 

function of calendar year, calendar month, and UT hour. 
While the trends are weak, at least some clear signatures 
are evident. First, the redshift completeness peaked in 
the early years of Swift (40% between 2005-2009) but 
has fallen since then (27% since 2010), probably because 
follow-up of “routine” bursts has become less common 
over the years and also because of smaller allocations of 
time to rapid-response spectroscopy. This trend is large 
enough that a year cut significantly improves the red¬ 
shift completeness of the sample; restricting ourselves to 
older bursts also has the benefit that the large amount 
of time elapsed has provided more opportunity for previ¬ 
ous host-galaxy campaigns to acquire ground-based data. 
We therefore restrict the base sample to events in the 
five-year 2005-2009 period. (Some post-2009 events are 
later included under a separate criterion; 1X71) 

A seasonal effect may also exist; the redshift com¬ 
pleteness drops significantly between January and March 
(26%) compared to the rest of the year (35%). We 
investigated whether this dependence may be different 
for Northern and Southern hemisphere bursts and found 
the trend to be present in both data sets. Despite the 
strength of this effect we elected not to cut on it, since 
the complementary benefits are less clear than for the 
other cuts employed. Nevertheless, considering the im¬ 
pact observed despite the crudeness of a raw seasonality 
cut, the trend we observed here provides good reason 
to consider more specific cuts, based on actual weather 
conditions at various sites, in the future. 

We see no clear secular dependence on the time of day 
and did not cut on this parameter. 

Figure [T] shows all GRBs through the end of 2014 in order to 
provide the most up-to-date view of the population, although our 
selection criteria were developed in late 2012. The trends are the 
same in either case. 
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Fig. 1.— Redshift recovery numbers and fractions for GRBs as a function of various observational parameters. We only consider 
rapidly-reported redshifts (in nearly all cases from absorption spectroscopy of the GRB afterglow) and do not include redshifts from host 
observations at late times. The blue histogram indicates all triggered Swift GRBs while the black histogram indicates Swift G RBs with 
known redshift. The lower panel shows the ratio of the two; error bars are calculated using the binomial Bayesian method of ICameronI 
II2011I ). Cuts employed by the SHOALS survey are indicated by thick black arrows. We also show the cuts employed by two other 
host-oriented uniform samples: TOUGH (green arrows) and BAT6 (orange arrows). 


3.2. Sun Angle 

Ideally, Swift’s pointing would always be centered close 
to 180 degrees away from the Sun, such that every burst 
detected by the BAT could be followed immediately by 
any telescope in the night hemisphere of Earth and re¬ 
main above the horizon until sunrise. In reality, opera¬ 
tional considerations (in particular the cooling require¬ 
ments of the XRT) require that the pointing direction of 
Swift is frequently far from this ideal and many GRBs are 
detected close in projection to the Sun, making ground- 
based follow-up challenging or impossible. 

In Figure [TJi-e we plot the redshift completeness 
against angular separation from the Sun as well as 
against separation in right-ascension only. It is clear that 
Sun angle has a large impact on redshift completeness, 
with about ^40% of bursts in the anti-Solar hemisphere 
successfully recovered but extremely poor completeness 
for bursts nearest the Sun (<10% within 60°). 

Previous surveys (TOUGH, BAT6) used a Sun dis¬ 
tance cut of > 55° and no right-ascension cut. How¬ 


ever, based on the plots above it is clear that somewhat 
more stringent cuts can significantly further improve the 
redshift completeness with only modest additional down- 
selection: the redshift efficiency remains very poor for 
events closer than ~ 75° or for bursts with a right as¬ 
cension difference of less than 5 hours. We employ a 
minimum right-ascension separation of 5 hours, and in 
principle an even more stringent cut could have been 
used at the expense of reducing the sample size. We 
make no formal cut on angular separation, but note that 
the right-ascension cut alone excludes all GRBs closer 
than 54° from the Sun. 

3.3. Lunation and Moon Angle 

While the location of the burst relative to the Sun has 
a large impact on the ability to follow up a burst, lu¬ 
nar considerations appear to be negligible. In Figure 
[TJ-g we plot the recovery rate against lunar illumination 
and moon distance, respectively. We see only weak, low- 
significance trends in either case, and so we do not use lu¬ 
nar information to restrict our sample. (While this result 
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may seem surprising, we note that bursts very close to 
the full Moon are also very nearly anti-Sun, and the im¬ 
proved observability may counteract the higher sky back¬ 
ground. We also note that bursts very close to the moon 
are automatically excluded from our sample, though not 
from FigurelU on the basis of the lack of XRT follow-up.) 

3 . 4 . Declination 

Bursts that occur close to the celestial equator can 
be observed by telescopes in both hemispheres, whereas 
those occurring closer to the poles are difficult or im¬ 
possible to follow from the opposing hemisphere and (in 
extreme cases) cannot be observed by equatorial-mount 
telescopes at all. On the other hand, high-declination 
bursts remain above the horizon for longer than equa¬ 
torial ones. In practice, the declination dependence of 
redshift completeness appears to be weak (Figure [T}r): 
no significant declination trend is seen in the Southern 
hemisphere, and while a fairly strong trend is evident 
in the Northern hemisphere it becomes much less sig¬ 
nificant after our preferred Ay and Tgo cuts ( il3.5Fil3.6p 
are applied, indicating that it may be coincidental. We 
therefore did not apply a declination cut. 

3 . 5 . Galaetic Latitude and Extinetion 

The position of a burst as projected against the plane 
of our own Galaxy strongly impacts our ability to recover 
its optical afterglow, mostly due to extinction by inter¬ 
stellar dust within the Milky Way (but crowding and 
source confusion with foreground stars may also play a 
role). This effect can be clearly seen in Figures [H and[T| 
where the recovery fraction is plotted against both the 
absolute Galactic latitude |&| a nd the foreground ex tinc- 
tion Ay (as calculated from the lSchlesrel et al.lll998l dust 
maps^"^). For extinction columns in excess of > 0.5 mag, 
and for latitudes less than |6| < 10 degrees, the recovery 
fraction plummets. 

We employ a maximum foreground extinction of Ay < 
0.5 mag (Schlegel) in SHOALS, which of course in prac¬ 
tice also cuts most events near the Galactic plane. This 
is the same criterion used by both TOUGH and BAT6. 

3 . 6 . Prompt Emission Properties: Fluence, Flux, and 
Duration 

The GRB prompt emission is not thought to depend 
directly on the properties of the circumprogenitor en¬ 
vironment, so incorporating prompt-emission properties 
in our selection should not bias the sample—although 
such selections will necessarily impact properties such as 
the intrinsic redshift distribution of the sample versus 
that of the Swift parent population (we expect brighter 
bursts to be closer, on average). More importantly, GRB 
prompt-emission brightness is also obser ved to be cor¬ 
related with GRB afterglow brightness ( Geh rels et all 
120081: iNvsewander et all 120091 iKarm et"^ 20101 1: partly 
because luminous bursts tend to release more kinetic en¬ 
ergy into their surroundings and produce brighter after¬ 
glows, and also simply because closer bursts (at fixed 
luminosity) will be brighter at all wavelengths. We there¬ 
fore expect that the brighter afterglows associated with 

We use the orig inal Schlegel et al. maps in preference to 
more recent work ('e.g.. lPeek fc Gravesll2nI(l : ISchlaflv fc Finkbein^ 
1201111 for consistency with earlier GRB host selection efforts. 


a population of bright GRBs should make their redshifts 
easier to measure. 

The impact of adding a prompt-emission cut on red- 
shift completeness wa s first clearly demonstrated by 
ISalvaterra et al.l (|2012h : the BAT6 sample reached ~80% 
redshift completeness even without any host-galaxy 
follow-up after application of a fairly stringent cut on 
the peak photon flux in addition to their observabil¬ 
ity criteria (after host follow-up this rose to ~90%). A 
flux/fluence cut is also desirable from the point of view of 
better understanding events with optical nondetections, 
since it is easier to establish l imits on the spectral index 
/3ox ll.lakobsson et al.l l2004 Ivan der Horst et all l2009ll 
and determine if the afterglow is obscured in the pres¬ 
ence of a brighter X-ray afterglow. Finally, imposition of 
an explicit cut on prompt emission properties ser ves to 
repla ce the complicated Swift triggering criteria (jBandl 
l2006ll with a simple well-defined threshold that can be 
more easily modeled for the purpose of measuring intrin¬ 
sic rates and redshift/luminosity distributions (|J6T]).^® 

We plot the dependence of the redshift completeness 
on the BAT 15-150 keV energy fluence and 1 s peak pho¬ 
ton flux in Figures HJc-l. The redshift completeness rises 
sharply for brighter bursts. Both fluence and peak flux 
show trends of similar magnitude, although in practice 
we found that using a fluence cut produced a slightly 
higher redshift completeness: after all other cuts were 
applied, the redshift completeness was 4% higher for a 
fluence cut than for a peak flux cut chosen to produce the 
same final sample size. A fluence cut is also expected to 
affect the redshift distribution less than a peak-flux cut 
would, since the cosmological dimunition in luminosity 
is offset by time dilation (see also HB.ll) . As a result we 
choose to cut on fluence (S'ls-isokeV > 10“® erg cm“^). 

We exclude short-duration bursts from the sample, re¬ 
quiring Tgo > 2 seconds. This was done not for its effects 
on redshift completeness, but rather to remove short- 
duration gamma-ray bursts from the sample, as these ap¬ 
pear to have a physically distinct origin associated with 
compact ob ject mergers and a very different host popu¬ 
lation fe.g.. iFong et al.l[2013l l. Strictly, this cut actually 
has no impact: our fluence cut strongly disfavors short 
events and removes all events with Tgo < 2 sec on its 
own. However, we also remove two events with Tgo > 2 
sec (GRBs 060614 and 080503) whose prompt-emission 
light curves resemble those of short-duration events with 
extended emission (while their exact origins are debated, 
these events have properties very different from ordi¬ 
nary long-duration GRBs: see e.g., iFvnbo et al.l 120061: 
Gai-Yam et al.1 120061 : iGehrels et al.l 120061: iPerlev et ahl 
2009bl for further discussion). 

Unlike the other properties discussed in this section, 
the BAT measurements are subject to measurement un¬ 
certainties and in principle values can change upon re¬ 
analysis. Our selection was conducted using values 
in the Swift GRB table, which is collated from val¬ 
ues reported in the BAT refined analysis GGN Gircu- 

We wish to emphasize that cutting on fluence may in principle 
affect the redshift distribution—but does so in a controlled way 
that improves, rather than harms, our ability to accurately measure 
the intrinsic burst rate. The BAT sample is always flux/fluence- 
limited; our cut replaces the complex on-board trigger criteria with 
a well-defined selection crite rion that we can correct for. This issue 
is discussed further in MM 
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lars (which in turn are based on analysis conducted 
by the BAT team 1-2 days after the event). Further 
analysis and improved measurements has been carried 
out o nly for a subset of these (e.g. ISakamoto et al.l 
unni). Since we wished to examine the entire Swift 
sample using the same data we used the GCN val¬ 
ues to establish our selection. Fortunately, the differ¬ 
ences are usually very small (a few percent) and the 
choice has negligible impact on the sample properties. 
Only four new GRBs would have entered the sample 
had we used the updated values o n bursts where the y 
are available: 060605 at z=3.78 (iFerrer o et al.l 12009(1 . 
060124 a t 2=2.300 (iFvnbo et al.l l2009bll . 090102 at 
z=1.547 (Ide Uearte Postigo et al.ll2009dl ^ . and 090728 
(unknown-z; possibl y an i?-dropout and at z ^ 6; 
iMelandri et alJ [2009( 1. No GRBs currently in the sam¬ 
ple would be dropped from it. 

3.7. Rapid Observability to Swift 

Some events which the Swift BAT detects cannot be 
rapidly observed by the XRT due to pointing consider¬ 
ations or technical difficulties. In such cases, the lack 
of a prompt X-ray position reduces some of Swift’s ad¬ 
vantages relative to previous satellites. Optical after¬ 
glow follow-up is more difficult because observers must 
identify an (often-faint) new source within the full 3' ra¬ 
dius initial BAT error circle, potentially causing after¬ 
glow identifications to be delayed or missed entirely and 
resulting in lower redshift completeness (a trend which 
can be seen in Figure [TJi). Furthermore, because later¬ 
time X-ray observations do not guarantee a detection, 
host-galaxy follow-up and redshift identification may be 
strictly impossible in some cases if early observations are 
not carried out. For both these reasons, we require a 
prompt XRT response (< 600 sec). 

3.8. Rapid Observability to Specific Ground Facilities 

Many of the factors described in earlier sections (Sun 
angle, etc.) are chosen based on their impact on ground- 
based observability. Of course, other factors we have not 
explicitly considered also affect whether or not a ground- 
based telescope can actually respond to a burst, includ¬ 
ing weather conditions, the instrument set(s) available 
that night, and whether or not the telescope was under 
maintenance. 

In principle it would be desirable to explicitly incor¬ 
porate site-dependent weather and telescope downtime 
at the most important facilities into our calculations as 
a formal sample cut. This may be a promising avenue 
for establishing future samples (in particular for more 
recent bursts, given the lower degree of redshift com¬ 
pleteness overall since 2010), but as this information is 
not readily available we did not consider it in designing 
our survey. However, we can be even more precise by 
considering whether or not several leading ground-based 
imaging response facilities actually were successfully trig¬ 
gered by Swift for early-time observations. Such an ap¬ 
proach has already been applied extensively by our team 
and others to produce uniform opt i cal aft e rglow samples 
(jCenko et al.l l20091j : iPerlev et al.l l2009d : iGreiner et al.l 

The change in status for GRB 090102 is not due to reanalysis 
but is the result of a typo in the Swift GRB table. 


I 2 OIIII and the same procedure can be applied to GRB 
hosts. 

Figure [TJr shows the recovery rate for events ob¬ 
served “rapidly” by various GRB-oriented follow-up fa¬ 
cilities: the Pa lomar 60-inch teles cope (P60; bursts 
from 2005-2012: iGenko et al.l[2009bl and work in prep.), 
the Gamma-Ray Bur st Optical/Near-inf rared Detector 
(GROND; 2007-2012: lGreiner et al.ll2011l and Kriihler et 
ah, priv. comm.), the Peters Automated I nfrared Tele¬ 
scope (PAIRITEL 2005-2010: lMorganll2014f) . the Robotic 
Optical Transient Search Experiment (ROTSE-III; 2005- 
2012; taken from a search of the GGN circulars), and the 
comb ination of Faulkes an d Liverpool telescopes (2005- 
2007: IMelandri et 11111200811 . compared to events observed 
by none of these facilities on a rapid timescale. Except 
in the case of the (small) Liverpool/Faulkes sample the 
triggering of early-time observations is, unsurprisingly, 
correlated with modest increases in recovery rate. 

Except possibly in the case of the ROTSE events 
(whose redshift completeness is lower than desired), these 
samples are not large enough for our purposes, even if 
applied with no other cuts. Therefore instead of requir¬ 
ing ground-based follow-up by one of these facilities, we 
treat rapid observations to either P60 or GROND as an 
alternative criterion allowing these bursts to bypass the 
three afterglow-observability requirements (year, hour- 
angle, and XRT-response time), such that a GRB nomi¬ 
nally failing one or more of these criteria can nevertheless 
enter the sample as long as either of these two telescopes 
observed it at early times. GROND or P60 observations, 
however, do not provide immunity to criteria not con¬ 
nected with ground observability, in particular the flu- 
ence and extinction cuts which remain in effect across 
the entire sample. 

To maximize o verlap with the TOUGH sample 
(|Hiorth et al.ll201^ . for which extensive host-galaxy data 
is already available, we offer a similar bypass to events 
within that sample (but, again, with the requirement 
that they must still satisfy the fluence cut.) 

3.9. Additional Host Observability Criteria 

Two additional criteria are added to ensure that host- 
galaxy follow-up is possible. First, we require that a 
< 2" (at 90% confidence) position be available, since 
uniquely identifying a host galaxy in an area larger than 
this carries significant risk of misidentification or ambi- 
gnity. Because all of our events were observed by the 
XRT at early times when the burst was bright (or by a 
sensitive ground-based optical facility) a position of this 
accuracy is almost always available, and this cut removes 
only a single event (GRB 0806I3B) from consideration. 
In principle this exclusion creates a very weak bias (if 
this event had an optical afterglow we would have in¬ 
cluded it: two other events, GRBs 07I02I and 061110B, 
have > 2" XRT position uncertainties but were included 
thanks to available optical/IR follow-up)—but given its 
small impact in practice we ignore this effect. 

We also require that the host position not be con¬ 
taminated by a bright foreground object in the form 
of a Galactic star or intervening galaxy. This decision 
is somewhat subjective (stars and galaxies contaminate 
many fields to various degrees, especially in the Spitzer 
imaging due to its large PSF), and in the case of a con¬ 
taminating galaxy requires some degree of deep follow- 
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up to have been obtained in the first place to recog¬ 
nize that the foreground galaxy is not the host. Events 
which we exclude on the basis of a contaminating bright 
star are GRBs 050716, 060923C, 071003, 080129, 080212, 
080229A, 080905B, and 101023A. Events excluded due 
to a contaminating galaxy are GRBs 080319C (for which 
ground-based spectroscopy and HST imaging shows the 
source closest to the afterglow location to be a superpo¬ 
sition of the host and a foreground system at z = 0.81) 
and 081028 (for which no ground-based imaging or spec¬ 
troscopy is available, but for which the observed IRAG 
magnitude of the source underlying the afterglow posi¬ 
tion is inconsistent with any galaxy at that redshift, so 
is probably one of th e foreground absorbers mentioned 
in iBerger et al.ll200§l . 

3.10. Summary and Sample Properties 

We summarize our final selection criteria^^ in Table [TJ 
Gharacteristics of the sample are presented in Table [H 
and a table of key properties of the sample relevant to our 
study (positions and redshifts; see subsequent sections) 
is presented in Table H) 

As these numbers indicate, our selection criteria were 
highly successful in isolating a large, well-observed sam¬ 
ple of Swift GRBs: 119 targets, 68% of which had 
redshifts measured before any late-time host follow-up 
had been conducted. Including redshifts from late-time 
follow-up (both preceding our efforts and including our 
host-galaxy campaign; see we have achieved a com¬ 
pleteness close to 90% (89% considering only secure spec¬ 
troscopic redshifts, 92% including photometric redshifts, 
94% including 77//-bracketed and a tentative single-line 
spectroscopic redshift, and 99% including upper limits) 
so far on a sample twice the size of previous efforts. 
Forthcoming observations will likely increase the com¬ 
pleteness even further. 

Approximately 20-30 events (^20% of the sample) 
can be classified as “dark”, with the exact number de¬ 
pending on the definition employed. 25 GRBs have 
no unambigous optical (0.3 — l.O^m) afterglow detec¬ 
tion reported in the GGN circulars or elsewhere, and 
in only four of these cases could this be readily at¬ 
tributed to lack of deep or early follow-up (see Table 
|4] for details.) On the other hand, 12 events with op¬ 
tical detections have red colors indicative of significant 
dust attenation, in some cases a great deal of dust at¬ 
tenuation (e.g.. Ay ~ 3 — 4 mag for GRBs 080607, 
090709A, and 100621A: iPerlev et ^ 120111 : iGenko et al.l 
I 2 OIOI: iGreiner et al.ll2013[l —these events probably lie in 
the same physical class as optically-undetected dark 
bursts, but optical detections were secured thanks to par¬ 
ticularly efficient follow-up and/or a very intrinsically lu¬ 
minous afterglow. 

Since the optical afterglow observations for many of 
the GRBs in our sample have not been thoroughly ana¬ 
lyzed beyond the quick reports given in the GGN circu¬ 
lars, it is of course possible that some additional events 
are (modestly) obscured without us being aware of it. 

While final for the purposes of this study, our sample is readily 
extendable by loosening or expanding the afterglow observability 
criteria. In hypothetical future cases where it becomes necessary 
to be specific, we will refer to the uniform sample established by 
the specific criteria outlined in Table[T]as the SHOALS09+ sample. 
In this paper we will simply refer to it as the SHOALS sample. 


Nevertheless, as our o bscured fraction is consistent with 
other recent estimates (fPerlev et al.ll2009(l: IGreiner et al.l 
I 2 OIII: ICovino et al.ll2013ll it is likely that the events we 
have identified constitute the large majority of all events 
that were dust-obscured in our sample. The afterglow 
properties of the sample (including quantitative metrics 
of darkness such as drrx: iJakobsson et al.l 1200^ will be 
revisited in more detail in forthcoming papers. 

4. OBSERVATIONS AND DATA ANALYSIS 

A primary goal of our survey is to produce high-quality, 
multi-filter SEDs for all host galaxies within the sam¬ 
ple, enabling the construction of rest-frame luminosity 
functions at any wavelength and the measurement of im¬ 
portant physical parameters (mass, SFR, etc.) via SED 
fitting. These efforts are still ongoing, and will be de¬ 
scribed in full in subsequent papers. Instead, we briefly 
outline our general observational strategy and its moti¬ 
vations. 

4.1. IRAC Observations 

Spitzer observations are key to our effort, extending 
the wavelength coverage by a factor of two relative to 
previous host galaxy surveys and providing access to 
physically distinct information to what is possible from 
ground-based observations alone. IRAG’s capabilities in 
its shortest-wavelength filters (3.6/im and 4.5/im) are 
undiminished even in its warm mission, and the in¬ 
strument remains sufficiently sensitive to detect typical 
galaxies out to z ^ 5. Furthermore, the luminosity of a 
galaxy at these wavelengths (which always probe wave¬ 
lengths redward of the Balmer break across this redshift 
range) is determined primarily by a single parameter 
(its stellar mass) with only modest dependence on age 
and extinction—providing a means of directly interpret¬ 
ing IRAG observations even without the supporting data 
we are amassing. (This contrasts with the situation at 
the rest-frame UV wavelengths probed by ground-based 
optical imaging, since a potentially very luminous galaxy 
can appear quite faint in these bands if it is heavily dust- 
obscured.) 

Observations of many targets in our sample were al¬ 
ready present in the Spitzer Legacy Archive; the remain¬ 
ing targets were observed as part of our Gycle 9 Large 
Program at 3.6 /rm with an exposure time chosen de¬ 
pending on the redshift (typically between 0.5-5 hours; 
see Paper II for details). Observations at 4.5 pm were 
also acquired for some targets, in particular for those at 
unknown redshift and for GRBs designated as “dark”, 
in order to provide better photometric redshift estimates 
and to model dust extinction in high-redshift galaxies. 

4.2. Optical and Near-Infrared Observations 

Spitzer observations alone provide an estimate of a 
galaxy’s total stellar mass, but do not constrain the other 
properties of a galaxy, including the nature of the young 
stellar population that (presumably) produced the GRB. 
In addition, because of Spitzer’s large PSF size (^2") 
it is not straightforward to uniquely identify the host 
galaxy based on Spitzer observations alone, even in pos¬ 
session of a precise afterglow localization: deep observa¬ 
tions approach the confusion limit, and it is not always 
clear whether an extended source at the afterglow loca- 
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TABLE 1 

Summary of Selection Criteria 


No. 

Type 

Short description 

AN’^ 

AC® 

1 

2 

GRB properties 

GRB properties 

Onboard Swift/BAT trigger before October 2012 w/XRT observations 
Tgo > 2s and not an SGRB+EE^^ 

-2 

+0.4% 

3 

GRB properties 

Sis-lsokeV > 10“® erg cm“2 

-69 

+9.6% 

4 

Host/afterglow visibility 

Galactic Ay < 0.5 mag 

-35 

+9.7% 

5 

Afterglow visibility 

Amenable to follow-up: (5a-i and 5a-ii and 5a-iii), OR 5b, OR 5c 

-143 

+17.1% 

5a-i 

Afterglow visibility 

XRT observations within 10 minutes 

-10 

+2.2% 

5a-ii 

Afterglow visibility 

Between 2005-2009 (inclusive) 

-36 

+8.1% 

5a-iii 

Afterglow visibility 

Sun hour angle separation >5 hours 

-13 

+5.7% 

5b 

Afterglow visibility 

Automatically triggered P60 within 1000 s or GROND within 1 hour 

+ IS 

+0.8% 

5c 

Afterglow visibility 

Satisfies TOUGH criteria H^ort^^^al^C^^') 

+11 

-1.4% 

6 

Host/afterglow visibility 

No known foreground star or galaxy contaminating the position 

-7 

+2.2% 

7 

Host visibility 

< 2" position available 

-1 

+0.6% 


^ Number of GRBs affected by this criterion, if applied after all other criteria; i.e., these numbers indicate the 
(negative of) the change in sample size if the criterion in question were removed and the sample reconstructed using 
all other criteria still in place. In the case of criteria 5b and 5c, numbers are positive since these criteria enable GRBs 
to be included despite failing one of the 5a criteria. 

Increase in redshift completeness after cutting on this criterion, if applied after all other criteria. This indicates the 
(negative of) our change in pre-host-followup redshift completeness if the given criterion was dropped. 

A short-duration GRB with extended emission. Only GRB 080503 and GRB 060614 fall under this category among 
targets not cut by other criteria. 


TABLE 2 

Sample Characteristics 


Total sample size 

119 

Number of “early” redshifts 

81 (68%) 

Number of redshifts to date 

110 (92%) 

Number with redshift limits 

118 (99%) 

Mean redshift 

2.18 

Redshift quartiles 

1.26, 2.06, 2.77 

Redshift range 

0.03-6.29 


tion represents an extended host-galaxy or a blend of the 
host galaxy and a foreground object. 

For both these reasons, we have obtained a large vol¬ 
ume of imaging at optical and near-lR wavelengths of all 
of our targets. Given the wide range in redshifts and lu¬ 
minosities of the galaxies targeted by our survey and its 
all-sky nature (requiring different observational facilities 
to cover the northern and southern regions, and observ¬ 
ing runs scattered throughout the year) this follow-up is 
necessarily heterogeneous and usually tailored to each in¬ 
dividual target. Typically, we try to obtain at least one 
deep (i?iim ~ 26) ground-based optical (rest-frame UV 
at z > 1) measurement and then obtain additional filters 
to the extent possible given the brightness of the target 
and the resources available; we also employ archival ob¬ 
servations from a variety of previous surveys and from 
the Gemini and VLT archives. 

Most observations were conducted using the Low R eso- 
lution Imaging Spectrometer (LRISi lOke et al.lll995D on 
Keck I, with substantial imaging also comi ng from the 
Gem ini Multi-Object Spectographs (GMOS; iHook et all 
l2004l) at Gemini-North and Gemini-South, the Optical 
System for Imaging and low-Intermediate-Resolution In¬ 
tegrated Spectroscopy (OSIRIS) at the Gran Telesco¬ 
pic Ganarias, the Inam ori Magellan Areal C amera and 
Spectrograph (IMACS; iDressler et al.l l2nilh at Magel¬ 
lan, the Focal Reducer and Low Dispersion Spectro¬ 
graph 2 (FORS-2) at the Very Large Telescope, the 
Wide Field Camera 3 (WFC3) aboard HST, and from 


the Gam ma-Ray Burst Optic al/Near-Infrared Detector 
(GROND; iGreiirer et al.l[200^ at the MPG 2.2m in La 
Silla. 

In general, we make use of standard reduction tech¬ 
niques and pipelines where available to produce a stacked 
image, then astrometrically align the stacked images 
against the Spitzer imaging. As the Spitzer PBCD 
imaging is by de fault astrometrically aligned against the 
2MASS catalog (ISkrutskie et al.ll2006D this effectively es¬ 
tablishes 2MASS as the astrometric reference system for 
the survey, and all positions reported in this work are 
therefore based on the 2MASS astrometric reference sys¬ 
tem. 

Observational efforts are still ongoing (although nearly 
complete for targets above 6 > —20°), so it is not yet pos¬ 
sible to present a complete catalog of imaging acquired by 
the survey: this will be presented in forthcoming work 
following the completion of this effort. At the present 
time, we have collected over 690 individual photometric 
data points on host galaxies within the sample, of which 
more than 510 represent detections (the remainder being 
upper limits; nonconstraining upper limits due to poor 
weather are excluded from these numbers). These are 
supplemented by additional photometry from the liter¬ 
ature. Every host galaxy except for one has at least 
one deep optical observation to supplement the IRAC 
data, and all except for 21 have at least one optical or 
infrared detection. Approximately half of the sample 
has numerous multicolor detections suitable for detailed 
characterization and modeling of the SED. We expect 
these statistics to improve modestly as our observational 
efforts wrap up during the coming year. 

4.3. Spectroscopic Observations 

A primary goal of the survey is to increase our spectro¬ 
scopic completeness as high as is possible to remove any 
bias associated with redshift measurement—in particu¬ 
lar that associated with dark bursts, but also potentially 
for bursts which do have detectable afterglows that are 
fainter than average and more difficult if not necessar¬ 
ily impossible to obtain absorption spectra of in time 
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(perhaps moderately extinguished bursts or those in low 
density media.) 

Many of the events in our sample without an ab¬ 
sorption redshift are prominent Swift dark bursts 
(or overlapped other uniformly-constructed surveys, 
such as TOUGH or BAT6) and as such a signifi¬ 
cant fraction of targets already had host-galaxy red- 
shifts i n the literature, primari ly from I Jakobsson et ^ 
(1201211: iSalyaterra et ahl (1201211: iKriihler et al.l ( ^1211 : 
iPerlev et al.l (|2013ll . or iKriihler et al.l ( 2015I1 . All re¬ 
maining sources showing bright (< 24 mag in any band) 
host-galaxy detections, and some fainter ones, were tar¬ 
geted for optical and/or NIR spectroscopy. For most 
of these observations we employ t he X-shooter spect ro- 
graph at the Very Large Telescope (jVernet et al.ll20lill . a 
medium resolution cross-dispersed echelle spectrograph 
simultaneously covering the wavelength range between 
0.3 fim and 2.5 /rm. For observations associated with 
our program (094.A-0593, PI S. Schulze), 4 exposures of 
900 s each were obtained in an ABBA nodding sequence, 
and reduced in a standar d manner using the X-shooter 
pipeline provided by ESO (|Goldoni et al.ll200^ and using 
our own routines. For No rthern-hemisphere ta rgets, we 
used LRIS or MOSFIRE (iMcLean et al.l[2?)lll on Keck 
I or NIRSPEG on Keck H and reduced the data using 
custom routines. 


4.4. Host Identification 

The number density of field galaxies on the sky is sig¬ 
nificant at the depths involved in our survey, so unam¬ 
biguously identifying the host galaxy requires as accu¬ 
rate and reliable a position of the originating GRB as 
possible. Wherever possible, we acquired the original 
target-of-opportunity imaging files showing the optical 
or near-infrared transient, aligned these against the late¬ 
time host imaging, and re-calculated the afterglow po¬ 
sition, which provides astrometric accuracy of typically 
0.4" or better (significantly smaller than the PSF of ei¬ 
ther the ground-based or Spitzer imaging). These images 
come from a number of sources, but common instruments 
are the imaging cam era on the Palomar 60-inch telescope 
(|Cenko et al.l 1200^ . the Nordic Optical Telescope, or 
publicly-available guider camera imaging from VLT or 
acquisition exposures from Gemini. In some cases we 
downloaded imaging cutout figures posted in the GCN 
circulars or in published papers, and re-calculate the po¬ 
sitions by aligning the cutouts to our host-galaxy imaging 
in a similar fashion. 

In cases where the original images are not available, 
we use astrometric coordinates published in the GCN 
circulars or in the literature, or supplied to us by others 
(in particular by D. Malesani and Y. Urata). We apply 
an astrometric offset (measured directly from the USNO 
and 2MASS catalogs) where necessary to translate from 
a USNO-aligned system to 2MASS: we assume published 
optical coordinates to be in a USNO-system if the astro- 
metric system is not stated explicitly (except for NIR 
imaging which we assume to be in the 2MASS system 
natively), although typically these offsets are quite small 
(^0.1-0.3") and do not dominate the uncertainty. Radio 
or millimeter coordinates are not offset. 

In some cases the only afterglow position available 
comes from Swift —typically from the XRT (as is often 
the case for dark GRBs), although in a few cases a Swift 


UVOT (UV-Optical Telescope: iRoming et al.l 120051) po¬ 
sition is available even though a ground-based position 
was not. Positions are taken from t he aut omated XRT 
analys es of iButler &: Kocevskil (|2007ll ^^ and lEvans et al.l 
(j2009t) ^^. which provide positions with a typical accu¬ 
racy of 1.5". UVOT positions are taken from the Swift 
GRB table. The coordinates are then shifted from their 
default USNO to the 2MASS frame, with the exception 
of positions noted as SDSS-aligned in the Butler tables 
which are not shifted. 

Since in the vast majority of cases we do have a sub- 
arcsecond position available, the probability of mistaken 
identification of the host galaxy of any individual well- 
localized GRB due to a chance foreground/background 
alignment is low—only approximately 1% of the sky is 
within 0.4" of a n unresolved i? < 26 mag galaxy (e.g., 
iHogg et ^Il997ll . Of course, in a large survey up to a few 
random alignments of this type would not be surprising, 
but would not significantly affect the results of the sur¬ 
vey. (Moreover, we can often identify and exclude them 
via a mismatch between host and afterglow redshifts: 
see 113.91) . The XRT-only, 1-2" positions are a source of 
somewhat greater concern, as the probability that a faint 
galaxy is present within a region of this size by chance 
is quite significant (^20%) and no prior redshift is avail¬ 
able. This is partially alleviated by the fact that XRT- 
only events are typically dark, and that dark GRBs tend 
to originate fro m hosts which are m uch brighter than av- 
erage (see e.g., iPerlev et al.ll2ni3l or IKriihler et al.ll201^ 
and have a low probability of chance association even 
considering the larger size of the XRT error circle. Only 
four sources have positional uncertainties greater than 
1" and host galaxies fainter than R = 25th magnitude 
(Fchance > 0.05): 050803, 050922B, 070621, and 070808. 
Furthermore, the first two of these sources are associated 
with blue-dropout (likely, z = 4 — 5) galaxies and the last 
is associated with an extremely red object—properties 
typical of optically-faint (dusty or high-z) GRB hosts 
but not common among galaxies selected randomly from 
the field. 


4.5. Host-Galaxy Photometry 

Magnitudes (or upper limits) for host galaxies observed 
as part of our survey are measured using aperture pho¬ 
tometry using a custom IDL wrapper around the aper 
photometric package included in the GSFC software li¬ 
brary, or are taken from the literature. Details of these 
procedures and the resulting photometry will be pre¬ 
sented in future papers (Paper H on the IRAC obser¬ 
vations is submitted concurrently with this paper), but 
are briefly summarized below. Photometry relevant to 
our photometric redshifts and upper limits is presented 
in Table p 

IRAC images at the depths relevant to our survey are 
at or near the confusion limit, and most of our host- 
galaxy targets show some degree of contamination from 
nearby foreground or background sources. To mitigate 
this contamination we employ an iterative, partially- 
automated PSF-fitting routine to model and subtract 
nearby objects identified from the Spitzer and ground- 
based imaging near the object and background apertures, 

http://butler.lab.asu.edu/Swift/xrt_pos.html 
http://www.swift.ac.uk/xrt_positions/index.php 
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Fig. 2. — New secure spectroscopic host-galaxy redshifts from our observations to date with VLT and Keck. Portions of the spectrum 
plotted in grey indicate wavelength regions affected by strong night-sky OH emission or by telluric absorption. (Many additional host 
spectroscopic redshifts are provided by our previously-published work and from the literature in general; see Table |4]) 


leaving an uncontaminated image of the host. The host 
magnitudes are then measured via aperture photometry, 
calibrated using the zeropoint data in the IRAC hand¬ 
book interpolated to the desired aperture using measure¬ 
ments of the instrumental PSF. 

Photometry of optical/NIR images is provided in a 
similar manner, minus the need (in nearly all cases) to 
subtract contaminating sources since the host galaxy is 
well-isolated. We calibrate relati ve to the latest release 
of the Sloan Digital Sky Survey (lAhn et al.ll2014fl wher¬ 
ever possible; otherwise, we provided our own secondary 
standards by observing the target on a photometric night 
alongside Landolt standards with a small telescope (P60 
for northern targets or GROND for southern targets). 
Magnitudes in non-SDSS filters are interpolated from the 
available griz magnitudes using the photometric trans¬ 
formation equations of Lupton (2005)^°. Calibration of 
HST photometry is performed using instrumental zero- 
points. 


5. NEW REDSHIFT CONSTRAINTS 
5.1. Host-Galaxy Redshift Measurements 


Eight new (not previously published^^) redshifts have 
so far been established by our program, either via the 
spectroscopic observations discussed in N4.3l or from pho¬ 
tometric redshift fitting to the photometry we have gath¬ 
ered so far. A summary of new redshifts provided by 
these methods is given below. Zoom-ins of the spec¬ 
tral lines are presented in Figure!^ and our photometric 
SEDs (showing the EaZy model for the best-fit redshift) 
are shown in Figure [2J the relevant data are presented in 
Table [51 

GRB050803 — The host galaxy is marginally (3-5 (t) 
detected in deep i?-band imaging from the VLT, in i- 
band imaging from the GTC, and with the F160W {H- 
band) filter on HST, but no other wavelength (although 
a hint of a marginal detection is evident in the Spitzer 
imaging). The lack of ( 7 -band detection is indicative of 
a dropout at a redshift of z ^ 4, although not an un¬ 
ambiguous one. Fitting the p hotometry using the E aZy 
photometric redshift software (jBrammer et al.ir2008f ) and 
disallowing passive low-redshift solutions indicates a red¬ 
shift of z ^ 4.3 (with ±0.4 redshift uncertainty at la 
confidence, although 1.9 < z < 4.9 is permitted at 95% 
confidence.) This redshift i s consistent with (b ut slightly 
higher than) the result of ISchulze et al.ll2015l using the 


Some of these redsh ifts are also reported in ISchulze et al.ll2015l 
and lKriihler et 50120151 . which were submitted close in time with 
http://www.sdss.org/drl2/algorithms/sdssUBVRITransform/#Lufltelfi200^ same investigators. 
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same data (but different photo-z software). 

GRB050922B — We observed this galaxy with LRIS {g 
and R filters) and with OSIRIS (i and z), complementing 
the existing R and Kg imaging from the TOUGH survey. 
The host galaxy is well-detected only in i, z, and in the 
Spitzer filters. A fit to the multi-filter SED (again ruling 
out passive low-z solutions) suggests that the redshift is 
high, 4.2 < z < 5.2 with a best-fit value o f z = 4.9, also 
in agr eement with the parallel analysis of iSchulze et al.l 

(Haii). 

GRB060204B — X-shooter observations of the host 
galaxy of this GRB reveal many strong (> 8a) emission 
lines corresponding to Ha, [O IIU 4959 ], [O IIUsoo?]: and 
[O IIA 3727 ], as well as probable H/3, at a common redshift 
of z = 2.3393. The velocity profile shows some structure 
with a prominent blue wing extending to about 300 km 
s“^. 

GRB061202 — X-shooter observations of the host 
galaxy show strong emission lines consistent with Ha and 
[O IIIa 5 oo 7 ] at a common redshift of z = 2.2543. Weaker 
emission is also seen at the locations of [O inA 4959 ] and 
H/3, and possibly [O IIA3727] (although the latter is in 
between two bright sky lines). 

GRB070223 — Spectroscopy of the host galaxy with 
LRIS shows two strong emission lines at wavelengths of 
A = 3198 A and A = 9801 A, corresponding to Lyman- 
a and the [O IIA3727] doublet at a common redshift of 
z = 1.6295. 

GRB080319A — We imaged this field extensively with 
Keck/LRIS (UBgRiZ) and Keck/MOSFIRE (YJHK). 
The host is well-detected in every band (and with IRAG). 
A fit to the SED indicates a photometric redshift of 
z = 2.43to.36. We also obtained iL-band spectroscopy 
of this target using MOSFIRE on the Keck I telescope 
on 2015-06-07 UT. A total of 20 exposures of 120 s each 
were obtained. A single, unresolved emission line is visi¬ 
ble in the subtracted and stacked 2D spectrum at a wave¬ 
length of 15153 A. This could be either [O IIIa 5007 ] at 
z = 2.0265 or Ha at z = 1.309, but the latter case is 
strongly ruled out by our photometric redshift. (An as¬ 
sociation with weaker lines such as H/3 or [O IIIA4959] is 
ruled out by the lack of additional line detections in the 
iJ-band spectroscopy.) We therefore infer a redshift of 
z = 2.0265. 

GRB080205 — A preliminary afterglow redshift of 
z ~ 4.0 was estimated by the UVOT team on the basis 
of an apparent drop out in the B filter in the ea rly-time 
UVOT photometry (|Oates fc Markwardtl l2008ll . How¬ 
ever, we strongly detect the host in our ground-based 
B—band and u-band imaging, ruling out a redshift this 
high. A photometric fit to our uBVRiz and Spitzer ob¬ 
servations imposes a maximum redshift of z < 3.08 with 
a best-fit redshift of z ^ 2.7 (with large uncertainties; 
2- — 2 71 + 0 - 25 \ 

— Z.l J-_Q ggP 

GRB 081210 — X-shooter observations of this tar¬ 
get show a strong emission line at 15336A and some 
weaker features. We identify the line as [O IIIA 5007 ] at 
z = 2.0631 on the basis of probable (4cr) detections of 
H/3 and [O IIIA4959] at a consistent redshift and the fact 
that if this were another strong line (e.g., Ha) other lines 
should be detected in clean regions of the spectrum but 
are not observed. At z = 2.0631 the wavelength cor¬ 
responding to Ha is in a region of moderately strong 



0.3 1 5 0.3 1 5 


Wavelength (gm) 

Fig. 3.— New photometric redshifts inferred by the survey. 
GRBs 050803 and 050922B show dropouts in the optical band 
indicative of a high redshift 2 : 4 — 5, although in the case of 

the former event this will need to be confirmed by deeper gi-band 
imaging. The remaining events show a weaker flux decrement in 
the C//B-bands indicative of the onset of the Ly-o forest and (for 
GRB 080319A) a Balmer jump in J-band at a consistent redshift. 
(The precise redshift of GRB 080319A was later fixed to 2 : = 2.0265 
via emission-line spectroscopy.) 

The model SEDs shown are from EaZy at the best-fit photometric 
redshift. 

telluric absorption. 

5.2. Redshift Upper Limits 

Nine sources have thus far eluded redshift measure¬ 
ment. However, even in most of these cases we can place 
upper limits due to the detection of optical afterglow, 
the detection of significant soft X-ra y absorption exces s 
above the foreground Galactic value ([GruDe et al.ll2007ll , 
or the detection of an optical host galaxy. 

GRB050128-We retrieved deep archival FORS-1 imag¬ 
ing of this source from the VLT archive; the data (associ¬ 
ated with program 075.A-0718(A)) were taken on 2005- 
05-11 UT and total 3080 seconds of exposure time in the 
i?-band filter. Within the XRT error circle we securely 
detect a host galaxy candidate in the combined stack. 
The detection in this band indicates z < 5.5. 

GRB050726-The host galaxy is not detected in a very 
deep VLT i?-band observation taken as part of the 
TOUGH project, so the host cannot be used to constrain 
the redshift. However, a detection of an early afterglow 
was reported in the initial UVOT U-band exposure, sug¬ 
gesting z < 3.5. Nondetections in subsequent UVOT 
B and [/-band observations suggest that the GRB may 
be a dropout and therefore close to this maximum red¬ 
shift, b ut it is also quite possible that the GRB faded 
rapidly (iPoole et al.l 1200511 — so only an upper limit can 
be placed. 

GRB070621-The (probable) host of this object, first 
identified in TOUGH, is very faint and detected only 
in B-band; while observations have been acquired in 
other filters with Keck and Gemini they do not reach 
the depths needed to detect the host. The afterglow was 
not detected at any optical band despite deep early-time 
imaging. The absorbing column inferred from the XRT 
spectrum is large, although not definitively so, and it 
places only a weak redshift limit of z < 5.5 (similar to 
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what is implied by the putative i?-band host detection). 

GRB070808- The afterglow of this event was detected 
only by XRT and the identity of the host galaxy itself 
is subject to some ambiguity given the coarse position; 
the nearest source to the center of either the DSS-refined 
or UVOT- refined XRT e rror c ircles was previously sug¬ 
gested by iHiorth et al.l (l2012ll as the most likely host 
candidate. This putative host galaxy is very red, well- 
detected in IRAC and in VLT ATg-band imaging but only 
marginally or not at all in the deep optical imaging ac¬ 
quired so far. Our optical photometry is not sufficient 
to establish a reliable photometric redshift on its own, 
but the red 3.6-4.5/rm color effectively rules out a low- 
redshift solution (requiring z > 0.5). This GRB also 
has one of the largest XRT Nh equivalent column ex¬ 
cesses measured in the entire sample, which in this case 
actually places a more constraining limit on the redshift 
than any properties of the host galaxy. Using the lower- 
limit Nh value measured from late-time XRT data by 
iBiitler fc Kocevskil ()2007ll of 4.06 x 10^^ cm ^ places a 
maximum redshift of z < 2.2 according to the empirical 
formula of lGrupe et al.l (|2007ll .^^ Treating the upper and 
lower limits together provides a crude redshift estimate 
of z = 1.35 ±0.85. 

GRB081128- A bright galaxy is located at the edge 
of the XRT error circle, although its centroid is signifi¬ 
cantly offset 1 ") from the optical afterglow position. 
A combination of multi-filter imaging and spectroscopy 
from both LRIS and X-shooter shows it to be an early- 
type galaxy at z = 0.27 with very little star-formation. 
In addition, a much fainter source is evident directly un¬ 
derlying the afterglow in LRIS g, and V band imaging 
(and marginally in i and z-band imaging with the same 
instrument, and with IRAC). While the nature of this 
system is not completely clear, the significant separation 
and color differential suggests that the bright galaxy is a 
foreground system unrelated to the GRB and the faint, 
blue object represents the true host galaxy. We fit the 
photometry of the fainter source using EaZy; while a con¬ 
sistent redshift of z = 0.32 is marginally favored higher 
redshift solutions (in particular, z ^ 3) are also credi¬ 
ble. We can place only an upper limit of z < 3.4 on the 
redshift. 

GRB100305A-No host galaxy is detected at the XRT 
position in our i maging (we note t hat the proposed opti¬ 
cal counterpart of lCucchiarall2010l is well outside the final 
XRT error circle and the source is still present in our own 
Keck imaging, so is not likely to be associated with the 
GRB). Two objects are just outside the edge of the error 
circle: one, to the southeast, is seen only in Spitzer and 
appears extended; the other, to the southwest, is also de¬ 
tected in our optical imaging. The uncertainty about the 
host identihcation and the lack of optical afterglow pre¬ 
cludes a definitive upper-limit based on these data. The 
excess X-ray column is also low, suggesting that this may 

To critically assess the IGrupe et al.l II2007I 1 relation we re¬ 
peated this procedure on the entire catalog of Swift events with 
spectroscopic redshifts, including dark bursts. With the exception 
of a handful of events with unrealistically large minimum Nh val¬ 
ues 10^® cm“^, probably due to intrinsic curvature or difficulty 
with the automated model-fitting; at face value these would ac¬ 
tually imply negative redshifts from the Grupe relation) we found 
only a single event with a redshift in excess of the predicted value: 
GRB 080207 (which exeeded it by only Az = 0.01). 


indeed be a high-redshift event that was “missed” on the 
basis of its faint afterglow. This is the only event in the 
sample without an upper limit on its redshift. 

GRB100802A-The onl y report ed detection of the op¬ 
tical afterglow is by P60 (|Genkoll201flf ). Only an r-band 
detection was reported originally. We stacked all of the 
imaging of the GRB acquired by the telescope that night 
and recovered detections in all four griz bands, and 
aligned these images against our late-time Keck imaging. 
A source is clearly detected in the Keck images at this 
location in R-band and (marginally) in R- and i-band, 
and with IRAG. The data are not sufficient to estimate 
a definite photometric redshift but place an upper limit 
of z < 3.1. 

GRB110709B-The host galaxy (first reported by 

1 Zanderer et al .112 0131) was observed with X-shooter (un¬ 
der VLT Program 090.A-0088), integrating for 8x900s 
using the AT-blocking filter. The reduced 2D frame does 
not show any strong lines over the spectral range, al¬ 
though a weak (5cr) line candidate is seen at an observed 
wavelength of 15479 A. The most likely identification of 
this feature is [O IIIasoot] at a redshift of z = 2.091, 
as other strong-line matches are ruled out ([O 11 ^ 3727 ] 
would be resolved, while Ha would imply emission lines 
of [O IIIa 5 oo 7 ] and [O 11 ^ 3727 ] within regions of the spec¬ 
trum where we have good sensitivity and do not ob¬ 
serve). We consider this assignment tentative pending 
future follow-up without the AT-blocking filter. Photo¬ 
metrically, the redshift can be limited only to z < 5.5 
(based on the HST F606W host detection). 

GRB120308A-The host is securely detected in our 
LRIS 5 -band and i?-band imaging from 2014-06-23 as 
well as in /-band LRIS imaging from 2014-05-27 and 
with IRAG. We attempted a photometic redshift fit to 
these data; this is unable to produce a lower limit on the 
redshift but does establish an upper limit of z < 3.7. 

6. RESULTS 

6 . 1 . Redshift Distribution 

The cumulative redshift distribution of our sample, as 
well as a binned fractional histogram, is plotted in Figure 
Remarkably, the redshift distribution of the sample 
is very similar to that of the overall Swift distribution. 
This suggests that, in spite of the many biases potentially 
affecting redshift measurement for a typical Swift GRB, 
the net impact on the redshift distribution is small. 

We investigated this more closely by calculating the 
redshift distributions of the dust-obscured GRBs and 
other GRBs separately (the most prominent instrinsic 
factor affecting redshift incompleteness is dust obscura¬ 
tion, which makes afterglow-based redshift determina¬ 
tion for ^20% of GRBs impossible in practice). The 
redshift distribution for obscured bursts is shown as a 

Here and in subsequent analysis we place photometric red¬ 
shifts (there are only six) at their best-fit values, and also use 
our best-guess redshift values for GRBs 070808 and 110709B. The 
remaining unknown-z bursts are omitted. Removing photomet¬ 
ric/insecure host redshifts would slightly decrease (by ~ 15%) the 
fraction of z ~ 3.5 — 5.5 events in the sample and have neglibile 
impact on the rest of the distribution; this would not influence our 
results/conclusions. If all unknown-z events were very close to thir 
maximum redshift, their inclusion would increase the fraction of 

2 ~ 3.5 — 5.5 events slightly (by ~ 20%) and likewise have neglible 
impact on the rest of the distribution ; again this would not affect 
our results or conclusions. 
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Fig. 4. — Redshift distribution of SHOALS GRBs (solid black 
line) compared to other GRB samples. Targets are mostly between 
0.6 < 5: < 4.0, with a few events at lower and higher redshifts; the 
overall distribution closely resembles that of known-redshift Swift 
GRBs as well as other complete GRB samples (e.g. TOUGH). 
Our modest fluence cut (Si5_i5okeV > 10~® erg cm~^) does not 
greatly affect the redshift distribution compared to that of Swift 
GRBs. 

red line in Figure SI Its distribution closely mirrors the 
general SHOALS population: a mild excess of GRBs is 
seen at 1.5 < z < 2.5 but a K-S test suggests this is 
not particularly significant (p = 0.15 comparing the ob¬ 
scured GRBs in the sample versus the remaining GRBs 
with redshifts.) This suggests that the fraction of cosmic 
star-formation that is obscured does not vary by a large 
amount with redshift, at least over the range in which 
we have reasonable number statistics (1 < z < 4). 

The redshift distribution is also very similar to that 
of the TOUGH sample (a similar unbiased sample—but 
with no fluence cut) indicating that, even allowing for the 
significant overlap between the samples, the fluence cut 
we employed has only a relatively minor impact. In con¬ 
trast, the BAT6 sample (which is cut fairly stringently 
on peak flux) shows a notable skew in its redshift dis¬ 
tribution towards lower redshifts. While our study is 
not contingent on the redshift distribution of our sample 
matching that of Swift GRBs overall (indeed, we would 
expect some differences based on the arguments in H3.6I) , 
this indicates that our final sample is nevertheless rea¬ 
sonably representative of the broader Swift population 
in redshift distribution. 

We observe a very small fraction of the sample at high 
redshifts. Only a single event is confirmed to be at 
z > 5.5 (GRB 050904 at z = 6.295), and among the 
9 GRBs in the sample with no measured redshift, all but 
one is limited to z < 5.5—so at most two events out of the 
119 in our sample can be at z > 5.5. This infrequency is 
qualitatively consist e nt with other comple t e studies (e.g., 
Perley et alJ 1200^ iFynbo et al.l l2009bl: [Greiner et alJ 
20111 : l.lakobsson et aLn2012l : ISalvaterra et al.l l2012t ) but 
even more constrai ning. Repeating th e Monte-Carlo 
analysis technique of iPerlev et al.l (l2009(J ) on our sample, 
we estimate that intrinsically between 0.3% and 5% of 
S > 10“® erg cm“^ Swift GRBs can originate at z > 5.5 
(95% confidence). 



Fig. 5.— Peak photon flux versus total fluence for Swift GRBs. 
Events in our sample are shown in yellow. Although the BAT 
triggers (approximately) according to photon count rate and our 
sample is cut in fluence, due to the correlation between the two 
parmeters our sample nevertheless is nearly complete in fluence: 
only the small triangular region to the right of the fluence cut 
and below the incompleteness threshold is expected to contain real 
bursts that would have passed our fluence cut if they had success¬ 
fully triggered the BAT. Diagonal lines are of constant F^/S and 
show the region containing nearly all long GRBs. 

While it is possible that we are preferentially missing 
GRBs at z > 6 in SHOALS relative to the all-Swift sam¬ 
ple due to the fluence cut, as discussed earlier in this 
section, the impact of a fluence cut at 10“® on the red¬ 
shift distribution appears to be small in practice. Our 
results therefore suggest that the high-z rate for Swift 
burs ts may be even lowe r than previously suspected (by 
e.g., iPerlev et alJ 1200931 : a few events per year. This 
underscores the challenges faced by recent efforts to ob¬ 
serve high-redshift GRBs, and helps to explain the re¬ 
cent paucity of confirmed high-z events. Likewise, it con¬ 
tributes to our finding that the intrinsic high-z GRB rate 
density is somewhat more modest than earlier estimates; 

m 

6.2. Sample Completeness With Respect to Prompt 
Emission Properties 

While we elected to cut on fluence^'* for the reasons 
detailed in 1 13.61 Swift's primary trigger mechanism is 
more strongly tied to a burst’s peak flux, and a difhcult- 
to-quantify incompleteness affects the satellite’s ability 
to detect and trigger on bursts with peak flux close to its 
threshold—complicating, in principle, any attempts to 
measure the intrinsic GRB rate. While a cut on peak flux 
at a value above the incompleteness level would largely 
eliminate this concern, the impact of our fluence cut is 
less straightforward: it is possible that Swift itself may 
have missed some GRBs whose fluence was above our 
cut level but whose peak flux was too low to trigger the 

Here and throughout the paper, we emphasize that the flux 
and fluence we refer to as measured refer only to the BAT 15-150 
keV band, not broad-band or bolometric values. Because GRBs 
occupy a wide range of intrinsic Speak there are many GRBs with 
bolometric fluences well above our cut that BAT is insensitive to, 
but because we consider BAT-band properties alone this does not 
impact our analysis or conclusions. 
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instrument (due to a particularly long and smooth light 
curve). 

Even so, we have reason to expect that the sample es¬ 
tablished by our fluence cut is nearly complete. In Figure 
El we plot the peak photon flux and total measured en¬ 
ergy fluence for all Swift GRBs (open circles) and for 
our SHOALS sample (filled circles).Unsurprisingly, 
these two parameters strongly correlate in a linear fash¬ 
ion but show scatter at the level of dex; the region of 
flux-fiuence space inhabited by most Swift long-duration 
GRBs (98%) is demarcated by the solid diagonals. The 
effective BAT triggering threshold is readily apparent in 
the data at approximately 0.4 ph cm“^ s“^. For bursts 
with very high fiuences, the BAT trigger sensitivity is 
clearly not an issue (e.g., no bursts with a fluence of 
S > 10“® erg cm“^ have a peak flux anywhere near the 
BAT incompleteness threshold: to a good approxima¬ 
tion, bursts this bright in the field of view should always 
trigger the telescope regardless of light curve shape, off- 
axis angle, etc.). For our chosen cut level of S' > 10“® 
erg cm“^, instrumental incompleteness is not completely 
negligible, and a few events which in principle may have 
been bright enough to meet our selection criteria could 
have been missed by the BAT. Based on the flux/fluence 
ratio distribution for brighter bursts where we are confi¬ 
dent the BAT sample is complete, we expect the number 
of such events is relatively small (^10 or less) and not 
likely to have a significant impact on the conclusions pre¬ 
sented in this paper. 

6.3. The Redshift-Dependent GRB Rate Density 

The observed GRB redshift distribution can be used 
to infer the intrinsic comoving GRB rate, provided that 
the sensitivity of the instrument and selection of the 
known-redshift sample are well-understood and the GRB 
luminosity function can also be inferred. As we argue 
above that the sample is intrinsically fluence-limited to 
a good approximation, these procedures can be applied 
to our sample also (using energetics-dependent quanti¬ 
ties in place of the more standard luminosity-dependent 
quantities). 

We first calculate the isotropic-equivalent energy re¬ 
leases for all GRBs in our sampl e, Uiso. 45 - 450 keV = 
5'i5-i50kev d-TT d\(l +z)~^k{z) (e.g.. iBlooin et al.ll200lD . 
Unlike most previous authors we do not make any at¬ 
tempt to estimate a bolometric or wide-bandwidth value, 
which is fraught with uncertainty considering that Swiff’s 
bandpass is narrow and does not usually contain the 
peak energy of the burst; the BAT-to-bolometric cor¬ 
rection can easily be orders of magnitude and it is prac¬ 
tically unconstrained by the Swift data. Instead, we cal¬ 
culate the energy release in the 45-450 keV rest-frame 
band only, corresponding to the window observed by the 
BAT (15-150 keV observer-frame) for an event at z ^ 2, 
close to the approximate sample redshift median and also 
where previously-reported divergences in the GRB-to- 
SFR ratio begin to become apparent. This makes the 
fc-correction much smaller and more reliable: specifically, 
its value is k{z) = [(I -I- z)l(l -I- 2)]^“^, where T is the 

Figure [5l and su bsequent figures/tables, show the updated 
IlSakamoto et al1l2011l'l values of flux and fluence for bursts where 
revised measurements of both values are available. Our results do 
not differ if the original GCN measurements are used. 



Fig. 6. — The comoving GRB rate density versus redshift com¬ 
pared to the comoving field-survey star-formation rate density (top 
panel), as inferred from our redshift and distribution (lower 
panel). Gurv es plot field-galaxy star-for ma tion histories from vari¬ 
ous sources (IM adau &: Dickins^ 12014 ] via IRobertson et al.ll2015l: 
Betnoozi_jtj)yj201j;JHo]3kins^£Beacom|J200fl;_ [Reddv &: Steidell 

200flnRobertson7rElli3l20T3 ~ l5ouwenr'SaLn2014ll . The scaling 

of the GRB rate (thick data points from the SHOALS sample, thin 
measurements from all S > 10~® erg cm“^ Swift GRBs; error bars 
denote 10-90% binomial confidence interval) is arbitrary, normal¬ 
ized against the star-formation rate at 2 ~ 2. The behavior of 
the GRB rate history is qualitatively similar to the star-formation 
rate history, but shows a modest (factor of ~ 5) excess at high 
redshifts (or equivalently, a depression at low redshift) compared 
to the most recent SFR density measurements. 

measured photon index over the BAT band (for bursts in 
which r is not well-constrained by the data we take the 
population median value of T = 1.5). This fc-correction 
is always small (exactly unity at z = 2 and within factor 
of 2 between I < z < 4 for any common value of T) and 
the uncertainty is even smaller (typically a few percent) 
so these measurements should be highly reliable. 

We assume an intrinsic isotropic-equivalent energetics 
rate function (f>{Eiso) following a single power-law and 
measure its power-law index by fitting the observed fiu- 
ence distribution of GRBs at moderate redshift (0.5 < 
z < 3.5) in the entire Swift sample above our fluence 
cut; we infer (j)Eiso oc We then count the ob¬ 

served number of GRBs (in both the SHOALS sample 
and for all S > I0“® Swift GRBs) with Eiso > 10^^ erg 
per redshift bin, and use the integral of the energetics 
function to correct the number in each bin for incom¬ 
pleteness in Uiso, scaling all bins to a common energy 
cutoff of Eiso > 10®^ erg. This number is then scaled in 
the usual way by dU/dz/(I -|- z) to convert the density 
in redshift to a comoving rate density. 

The resulting redshift-dependent rate (top panel of 
Figure Eland Table [3]) shows a broad peak atz~ 1.5 —3.5 
and a modest decline towards lower and higher red- 
shifts. This behavior is consistent with previous studies 
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TABLE 3 

GRB Rate Measurements 


•^min 

^max 

FJ - ' 

-^iso,sens 

SHOALS 

Gcomplete® PGRB,51® 

All Swift 

N'’ PGRB,51® 

0.10 

0.50 

0.90 

0.91 

2 

0 

5 

0 

^•^o_ 0 .io 

0.50 

1.00 

3.19 

2.83 

13 


26 

-‘^•^°-0.29 

1.00 

1.50 

6.22 

5.17 

11 

1.80t°;«i 

28 

Q r)Q”t"0.58 
^•^°-0.44 

1.50 

2.00 

9.60 

7.65 

16 

3-95tJ:®® 

31 

001 +0.91 

2.00 

2.65 

14.19 

10.8 

27 

8-34ll:li 

45 


2.65 

3.50 

20.23 

14.9 

17 

6.90tlil 

38 

'•'^'-1.28 

3.50 

4.50 

27.19 

19.5 

7 


15 

4 /)4+1-73 

4.50 

6.00 

37.15 

25.8 

5 

0 '7C+2.98 
'J- <^-1.39 

9 

o.ou_i Q3 

6.00 

8.00 

49.54 

33.5 

1 

-*^•-*-^-0.53 

2 


8.00 

10.00 

61.03 

40.4 

0 

u.uu_o 00 

0 

u.uu_o 00 


^ Eiao sensitivity threshold for inclusion of a GRB in the count, in units 
of 10®^ erg as measured in the 45-450 keV rest-frame band. 

Completeness correction applied to scaling observed GRB counts to a 
common threshold of Eiso > 10®*^ erg. 

Number of GRBs observed with Eiso > 7?iao,sens 

On-axis luminous (i?iso, 45-450 > 10®^ erg) GRB rate density over this 
redshift interval, in units of 10“® yr~^ Mpc“^. 


of this type (e.g.. iKistler et al, 
2010at fWanderman fc Pirad 


2012t [Jakobsson et al.l l20Ta 


[2rM llin^ iButler et al 


2010t iRobertsqn fc Ellii 


Salvaterra et al.l I2012D — 


but, unlike these previous studies (wit h the exceptions 
for th e si gnificantly smaller sam ples of I Jakobsson et ^ 
I2OI2I and ISalvaterra et al.ll201^ our results are not lim¬ 
ited by systematics associated with the highly incomplete 
redshift distribution. 

It is somewhat surprising that the complete, unbiased 
redshift distribution produced from our work is so similar 
to the distribution inferred from previous studies based 
(largely) on samples drawn from afterglow redshifts in 
the literature, which favor low-extinction sightlines. As 
we mentioned earlier ( J6.3I) this suggests a relatively uni¬ 
form fraction of obscured star-formation across most of 
cosmic history (our sample provides good number statis¬ 
tics between approximately 0.5 < z < 4). It also in¬ 
dicates that high-z GRBs, despite being uncommon (in¬ 
trinsically and observationally) and difficult to follow-up, 
are being identified with comparable efficiency as their 
low-redshift counterparts. 

7. A REDSHIFT-DEPENDENT GRB EFFIGIENCY? 

Since GRBs originate from massive stars, the com¬ 
parison of the GRB rate history and the star-formation 
rate history as derived by traditional galaxy sur¬ 
vey methods (e.g.. iMada u et al.lll998t iReddv fc St eidell 


20091 iHopkins fc Beacom l2006t iBehroozi et al.l l2013t 


Bouwens et al. I I20l4 see lMadau fc DickinsonI l2014l for 

a review) imposes a constraint on the degree to which 
the GRB rate is affected by other factors, as well as 
an independent check on the star-formation rate his¬ 
tory itself. Many previous studies have quantitatively 
compared these two relations in detail; most of these 
(iDaigne et al.ll200'^ iLe fc Dermedl2007t iGuetta fc PiranI 


aig- .. - -- .. . - 

200^ l^lvaterra fc Ghincarinill2007t lYiiksel et al. l2008t 

Kistler et al.ll2008l: iVirgili et al.ll201li iRobertson fc Ellis! 


2012t iLien et al.l 12014 c.f. Wanderman fc PirarJ l2010t 


Elliott et al.l 2012r Yu et al.l 20151) have concluded that 

the GRB rate at high redshift is significantly higher than 


what would be inferred from the galaxy survey-inferred 
star-formation history. This suggested either that the 
field surveys were insufficiently accounting for the num¬ 
ber density of low-luminosity galaxies they do not de¬ 
tect, or that the cosmic GRB production efficiency e(z) = 
Pgrb(-z)/psfr(- 2) is not constant (due to e.g. metallicity 
enrichm ent suppressing the GRB rate at lowe r redsh ift): 
see e.g.. iJakobsson et al.l (I2012D ; iTrenti et all (l2013l) . 

We carried out this exercise with our own observations 
as well, comparing our (selection-unbiased) rate distri¬ 
bution to the most recent estimate of the star-formation 
rate density out to very high re dshifts using the re¬ 
cent g alaxy luminosity functions of IMadau &: DickinsonI 
(I2014I). integrat e d dow n to L = 10“^L* galaxies by 
IRobertson et al.l (120151) . We confirm the high-redshift 
excess (or, equivalently, a low-redshift deficiency) in the 
GRB rate relative to the UV-inferred star-formation rate: 
normalizing^® the two curves at z = 2, the z ^ 5 GRB 
rate is in excess of the star-formation rate by a factor of 
2-3 and the z < 0.5 GRB rate is below it by a factor of 
2-5. 

This provides further support to the notion that the 
cosmic GRB efficiency may vary over time. We empha¬ 
size, however, that while the deviation is significant, it is 
also relatively modest in magnitude: a factor of ^5-10 
across the entire span of cosmic history from z = 5 to 
z ^ 0. In fact, perhaps the more salient conclusion to 
be drawn from Figure [6] is that strong evolution in the 
GRB-to-SFR ratio is not observed. This argues that the 
GRB rate’s dependency on metallicity must have only 
low-to-moderate impact on the cosmological rate, in dis¬ 
agreement with models i mplying strong variation s (e.g., 
the single-star models of iLanger fc Normanll2006l which 
require approximately Z < 0.1 Zq, would imply a varia¬ 
tion by more than a factor of 40 over this period). Ear- 

Note that the relative normalization of the two curves is ef¬ 
fectively arbitrary since we do not know the fraction of stars which 
explode as GRBs or the beaming correction to better than an order 
of magnitude. 
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Her analyses using: non-uniform afterglow-based samples 

(iRobertson fc Eliyl2012l : [Trenti et ^ 20131 : iHao fc Ynanl 

l2013f) have reached similar conclusions but were limited 
by systematics associated with incomplete redshift mea¬ 
surement; our work conhrms that these conclusions hold 
within an unbiased sample. 

Our result is also in agreement with independent inves¬ 
tigations of the variations of the GRB rate based on the 
properties of t he hosts themselves. Recent emission-line 
studi es (e.g., iGraham fc Fruchteil 120131 : iKriihler et al.l 
l2015ll show that GRBs form readily at moderate metallic- 
ities (~0.5Zq or more) that are characteristic of typical 
star-forming galaxies throughout most of the Universe’s 
history. And the host stellar mass distribution we in¬ 
fer from IRAG observations of the host galaxies in our 
sample (Paper 11) similarly suggests that a metallicity 
threshold for GRB production is likely present but is rel¬ 
atively high, ^ Zq. 

8. SUMMARY 

We have dehned a new legacy sample of gamma-ray 
bursts and host galaxies spanning nearly all of cosmic 
history: 119 Swift GRBs at 0.03 < z < 6.29 drawn 
from the Swift GRB catalog using a series of observabil¬ 
ity cuts, plus a requirement that the BAT fluence exceed 
S > 10“® erg cm“^. Gombining publicly-available af¬ 
terglow and host-galaxy redshift measurements with our 
own host-galaxy campaign, we present redshifts for 110 
(92%) of these bursts, or 112 (94%) if we include the 
tentative redshift of 110709B and the lower-|-upper lim¬ 
its on 070808. This is by far the largest, and the most 
redshift-complete, sample of its type, and provides the 
most up-to-date and unbiased view of cosmic history as 
seen by GRBs. Only one event lacks a redshift upper 
limit and only one event with measured redshift is at 
z > 5.5. 

Mapping our redshift distribution to a comoving rate 
density to measure the evolution of the GRB rate with 
cosmic time, we measure a rise in the GRB rate from 
z > 6 to z 2, followed by a drop of an order of magni¬ 
tude from z ~ 2 to the present time—the same pattern 
seen by traditional metrics of the cosmic star-formation 
rate density. Quantitatively comparing the GRB rate 
history to the cosmic star-formation history, we hnd a 
modest excess in the GRB rate (versus SFRD) at high 
redshift compared to lower redshifts. Consistent with 
previous work, this suggests that the cosmic GRB effi¬ 
ciency was higher in the first few billion years of cosmic 
history relative to today and provides support to the 
idea of a metal-dependent progenitor, but the modest 
degree of this variation rules out models requiring an ex¬ 
clusively very-metal-poor (e.g., < O.IZq) environment. 
The small number of high-redshift GRBs in the sample 
places strong limits on the fraction of high-z bursts de¬ 
tected by Swift and on the intrinsic GRB rate at high 
redshifts. 

Addressing the GRB rate-evolution question in 
detail—and actually applying our GRB population to 
address broader questions in astronomy—requires more 
than just redshifts. In particular, examination of the af¬ 
terglows and (especially) host galaxies of these events is 
needed to study the galaxy population giving rise to the 
GRBs in our sample at each redshift, providing an in¬ 
dependent test of factors controlling the GRB rate and 


a means to explore directly the importance and nature 
of the dusty, low-luminosity, and high-redshift popula¬ 
tions uniquely probed by GRBs. We are collecting and 
analyzing these observations under the programs intro¬ 
duced here, and the direct study of the hosts within our 
sample will serve as the subject of all remaining papers. 
As the first large, thorough, highly-complete, and multi¬ 
band survey of an unbiased GRB host sample, SHOALS 
will enable unprecedented insight into the GRB rate and 
progenitor as well as a unique perspective into galaxy 
evolution and cosmic history. 
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TABLE 4 

The shoals GRB sample 


GRB 

Seat® 


dec^ 

unc.^ 

type.'i 

OT?e 

dusty?^ 

ref.' 

z 

early?s 

origin^ 

ref.' 

050128 

51 

14:38:17.68 

-34:45:55.4 

1.70 

X 



1 

<5.5 

N 

host 

* 

050315 

32 

20:25:54.17 

-42:36:02.2 

0.28 

o 

Y 


2 

1.9500 

Y 

AG 

3 

050318 

11 

03:18:51.01 

-46:23:44.0 

0.60 

o 

Y 


4 

1.4436 

Y 

AG 

3 

050319 

13 

10:16:47.94 

+43:32:53.5 

0.40 

o 

Y 


5 

3.2425 

Y 

AG 

6 

050401 

83 

16:31:28.80 

+02:11:14.0 

0.43 

o 

Y 

Y 

2 

2.8983 

Y 

AG 

6 

050525A 

151 

18:32:32.59 

+26:20:22.3 

0.31 

o 

Y 


2 

0.606 

Y 

host 

7 

050726 

20 

13:20:11.95 

-32:03:51.2 

0.32 

u 

Y 


2 

<3.5 

N 

AG 

* 

050730 

24 

14:08:17.10 

-03:46:17.6 

0.37 

o 

Y 


8 

3.9693 

Y 

AG 

6 

050802 

22 

14:37:05.84 

+27:47:12.4 

0.60 

o 

Y 


* 

1.7102 

Y 

AG 

6 

050803 

22 

23:22:37.86 

+05:47:08.0 

1.40 

X 

N 

Y 

1 

4 3+0.6 
‘i.o 2.4 

N 

host 

* 

050814 

20 

17:36:45.39 

+46:20:21.2 

0.30 

o 

Y 


9 

5.3 

Y 

AG 

10 

050820A 

34 

22:29:38.10 

+19:33:36.8 

0.31 

o 

Y 


8 

2.6147 

Y 

AG 

6 

050822 

25 

03:24:27.22 

-46:02:00.1 

0.70 

X 



11 

1.434 

N 

host 

12 

050904 

52 

00:54:50.88 

+14:05:09.5 

0.40 

o 

Y 


* 

6.295 

Y 

AG 

13 

050922B 

24 

00:23:13.38 

-05:36:17.5 

1.70 

X 

N 


1 

4 Q+O.^ 

N 

host 

* 

050922C 

16 

21:09:33.08 

-08:45:30.3 

0.30 

o 

Y 


14 

2.1995 

Y 

AG 

6 

051001 

18 

23:23:48.72 

-31:31:23.6 

1.50 

X 

N 

Y 

11 

2.4296 

N 

host 

15 

051006 

12 

07:23:14.10 

+09:30:19.4 

1.50 

X 

Y 

Y 

1 

1.059 

N 

host 

12 

060115 

17 

03:36:08.32 

+17:20:42.8 

0.31 

o 

Y 


8 

3.5328 

Y 

AG 

6 

060202 

22 

02:23:22.94 

+38:23:03.9 

0.50 

I 

N 

Y 

* 

0.785 

N 

host 

16 

060204B 

29 

14:07:14.89 

+27:40:36.2 

1.10 

o 

Y 


17 

2.3393 

N 

host 

* 

060210 

76 

03:50:57.38 

+27:01:34.2 

0.60 

o 

Y 

Y 

18 

3.9122 

Y 

AG 

6 

060218 

65 

03:21:39.68 

+16:52:01.9 

0.28 

o 

Y 


8 

0.0331 

Y 

host 

19 

060306 

22 

02:44:22.92 

-02:08:54.1 

1.30 

X 

N 

Y 

11 

1.559 

N 

host 

20 

060502A 

23 

16:03:42.62 

+66:36:03.0 

0.40 

o 

Y 


* 

1.5026 

Y 

AG 

6 

060510B 

40 

15:56:29.48 

+78:34:12.1 

0.20 

o 

Y 


* 

4.9 

Y 

AG 

21 

060522 

11 

21:31:44.84 

+02:53:09.7 

0.42 

o 

Y 


8 

5.11 

Y 

AG 

22 

060526 

12 

15:31:18.34 

+00:17:04.9 

0.16 

o 

Y 


8 

3.2213 

Y 

AG 

6 

060607A 

26 

21:58:50.40 

-22:29:47.1 

0.38 

I 

Y 


8 

3.0749 

Y 

AG 

6 

060707 

16 

23:48:19.06 

-17:54:17.3 

0.37 

o 

Y 


8 

3.4240 

Y 

AG 

6 

060714 

29 

15:11:26.41 

-06:33:58.3 

0.43 

o 

Y 


8 

2.7108 

Y 

AG 

6 

060719 

15 

01:13:43.71 

-48:22:51.0 

0.29 

o 

Y 

Y 

8 

1.5320 

N 

host 

15 

060729 

26 

06:21:31.80 

-62:22:12.3 

0.27 

o 

Y 


8 

0.5428 

Y 

AG 

6 

060814 

148 

14:45:21.31 

+20:35:10.5 

0.18 

o 

IR 

Y 

8 

1.9229 

Y 

host 

15 

060908 

28 

02:07:18.41 

+00:20:31.3 

0.42 

o 

Y 


8 

1.8836 

Y 

AG 

6 

060912A 

14 

00:21:08.14 

+20:58:17.4 

0.35 

o 

Y 


8 

0.937 

Y 

host 

23 

060927 

11 

21:58:12.01 

+05:21:48.6 

0.19 

o 

Y 


8 

5.467 

Y 

AG 

24 

061007 

450 

03:05:19.58 

-50:30:02.3 

0.33 

o 

Y 


8 

1.2622 

Y 

AG 

6 

061021 

30 

09:40:36.12 

-21:57:04.8 

0.30 

o 

Y 


8 

0.3463 

N 

AG 

6 

061110A 

11 

22:25:09.84 

-02:15:31.4 

0.37 

o 

Y 


8 

0.7578 

Y 

AG 

6 

061110B 

14 

21:35:40.39 

+06:52:34.0 

0.24 

o 

Y 


8 

3.4344 

Y 

AG 

6 

061121 

139 

09:48:54.55 

-13:11:42.9 

0.36 

o 

Y 


8 

1.3145 

Y 

AG 

6 

061202 

35 

07:02:06.09 

-74:41:54.7 

1.40 

X 



1 

2.253 

N 

host 

* 

061222A 

81 

23:53:03.41 

+46:31:58.6 

0.30 

I 

IR 

Y 

25 

2.088 

N 

host 

26 

070110 

16 

00:03:39.27 

-52:58:27.2 

0.30 

o 

Y 


8 

2.3521 

Y 

AG 

6 

070129 

30 

02:28:00.94 

+ 11:41:04.1 

0.31 

o 

Y 


8 

2.3384 

N 

host 

15 

070223 

19 

10:13:48.39 

+43:08:00.7 

0.40 

I 

Y 


27 

1.6295 

N 

host 

* 

070306 

55 

09:52:23.30 

+10:28:55.2 

0.31 

o 

IR 

Y 

8 

1.4959 

Y 

host 

28 

070318 

26 

03:13:56.81 

-42:56:46.1 

0.35 

o 

Y 


8 

0.840 

Y 

AG 

29 

070328 

91 

04:20:27.73 

-34:04:00.5 

1.40 

X 



1 

2.0627 

N 

host 

20 

070419B 

75 

21:02:49.77 

-31:15:49.0 

0.39 

o 

Y 


8 

1.9588 

N 

host 

15 

070508 

201 

20:51:11.70 

-78:23:05.1 

0.40 

o 

Y 

Y 

* 

0.82 

Y 

host 

30 

070521 

81 

16:10:38.61 

+30:15:21.9 

1.40 

X 

N 

Y 

1 

2.0865 

N 

host 

20 

070621 

44 

21:35:10.09 

-24:49:03.1 

1.40 

X 

N 

Y 

1 

<5.5 

N 

host 

* 

070721B 

36 

02:12:32.96 

-02:11:40.8 

0.38 

o 

Y 


8 

3.6298 

Y 

AG 

6 

070808 

13 

00:27:03.36 

+01:10:33.9 

1.50 

X 

N 

Y 

1 

1.35+0.85 

N 

A^h.x 

* 

071020 

23 

07:58:39.78 

+32:51:40.4 

0.35 

I 

Y 


31 

2.1462 

Y 

AG 

6 

071021 

14 

22:42:34.30 

+23:43:06.2 

0.60 

o 

IR 

Y 

32 

2.4520 

Y 

host 

15 

071025 

73 

23:40:17.07 

+31:46:42.8 

0.35 

I 

Y 

Y 

33 

6 c 
+ 1 
00 

N 

AG 

34 

071112C 

30 

02:36:50.95 

+28:22:16.8 

0.41 

o 

Y 


5 

0.8227 

Y 

AG 

6 

080205 

20 

06:33:00.63 

+62:47:31.7 

0.50 

u 

Y 


35 

9 79+9-^4 
0.74 

N 

host 

* 

080207 

61 

13:50:02.97 

+07:30:07.3 

0.50 

X 

N 

Y 

36 

2.0858 

Y 

host 

15 

080210 

18 

16:45:04.01 

+13:49:35.6 

0.60 

u 

Y 


35 

2.6419 

Y 

AG 

6 

080310 

23 

14:40:13.80 

-00:10:30.7 

0.40 

o 

Y 


37 

2.4274 

Y 

AG 

6 

080319A 

44 

13:45:20.01 

+44:04:48.4 

0.70 

o 

Y 


38 

2.0265 

N 

host 

* 

080319B 

850 

14:31:40.99 

+36:18:08.7 

0.30 

o 

Y 


* 

0.9382 

Y 

AG 

6 

080325 

49 

18:31:34.23 

+36:31:24.8 

0.30 

I 

IR 

Y 

39 

1.78 

N 

host 

40 

080411 

265 

02:31:55.21 

-71:18:07.3 

0.50 

u 

Y 


35 

1.0301 

Y 

AG 

6 

080413A 

35 

19:09:11.75 

-27:40:40.4 

0.40 

o 

Y 


* 

2.4330 

Y 

AG 

6 

080413B 

33 

21:44:34.66 

-19:58:52.4 

0.50 

u 

Y 


35 

1.1014 

Y 

AG 

6 

080430 

12 

11:01:14.76 

+51:41:08.0 

0.40 

o 

Y 


41 

0.767 

Y 

AG 

42 

080603B 

25 

11:46:07.67 

+68:03:39.8 

0.30 

u 

Y 


35 

2.6892 

Y 

AG 

6 

080605 

134 

17:28:30.04 

+04:00:56.0 

0.30 

o 

Y 

Y 

* 

1.6403 

Y 

AG 

6 

080607 

247 

12:59:47.21 

+15:55:10.5 

0.40 

o 

Y 

Y 

* 

3.0368 

Y 

AG 

6 
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TABLE 4 

The shoals GRB sample 


080710 

14 

00:33:05.63 

-1-19:30:05.4 

0.30 

0 

Y 


* 

0.8454 

Y 

AG 

6 

080721 

141 

14:57:55.84 

-11:43:24.6 

0.40 

0 

Y 


5 

2.5914 

Y 

AG 

6 

080804 

38 

21:54:40.18 

-53:11:04.8 

0.35 

0 

Y 


* 

2.2045 

Y 

AG 

6 

080805 

27 

20:56:53.45 

-62:26:40.0 

0.60 

0 

Y 


43 

1.5042 

Y 

AG 

6 

080810 

46 

23:47:10.49 

-1-00:19:11.3 

0.30 

0 

Y 


* 

3.3604 

Y 

AG 

6 

080916A 

42 

22:25:06.23 

-57:01:22.9 

0.40 

0 

Y 


* 

0.6887 

Y 

AG 

6 

080928 

24 

06:20:16.83 

-55:11:58.7 

0.40 

0 

Y 


* 

1.6919 

Y 

AG 

6 

081008 

43 

18:39:49.86 

-57:25:53.0 

0.60 

U 

Y 


35 

1.967 

Y 

AG 

44 

081029 

21 

23:07:05.35 

-68:09:19.7 

0.50 

u 

Y 


35 

3.8479 

Y 

AG 

45 

081109A 

40 

22:03:09.61 

-54:42:40.1 

0.35 

0 

Y 

Y 

46 

0.9787 

N 

host 

15 

081118 

12 

05:30:22.18 

-43:18:05.1 

0.60 

0 

Y 


47 

2.58 

Y 

AG 

48 

081121 

42 

05:57:06.17 

-60:36:09.8 

0.60 

u 

Y 


35 

2.512 

Y 

AG 

49 

081128 

23 

01:23:13.10 

-1-38:07:38.7 

0.20 

0 

Y 


* 

<3.4 

N 

host 

* 

081210 

19 

04:41:56.20 

-11:15:26.8 

0.64 

u 

Y 


35 

2.0631 

N 

host 

20 

081221 

189 

01:03:10.16 

-24:32:51.6 

0.25 

I 

IR 

Y 

* 

2.26 

N 

host 

50 

081222 

52 

01:30:57.60 

-34:05:41.6 

0.50 

u 

Y 


35 

2.77 

Y 

AG 

51 

090313 

15 

13:13:36.20 

-1-08:05:49.6 

0.70 

0 

Y 


52 

3.375 

Y 

AG 

53 

090404 

31 

15:56:57.52 

-1-35:30:57.5 

0.40 

R 

N 

Y 

54 


N 

host 

16 

090417B 

23 

13:58:46.58 

-1-47:01:04.8 

1.00 

X 

N 

Y 

11 

0.345 

N 

host 

55 

090418A 

47 

17:57:15.16 

-1-33:24:20.9 

0.40 

0 

Y 


56 

1.608 

Y 

AG 

57 

090424 

218 

12:38:05.08 

-1-16:50:15.0 

0.60 

0 

Y 


* 

0.544 

Y 

host 

58 

090516A 

90 

09:13:02.60 

-11:51:15.0 

0.40 

0 

Y 


* 

4.109 

Y 

AG 

59 

090519 

12 

09:29:07.00 

-1-00:10:48.9 

0.60 

0 

Y 


60 

3.85 

Y 

AG 

61 

090530 

11 

11:57:40.49 

-1-26:35:37.7 

0.40 

0 

Y 


* 

1.266 

N 

host 

62 

090618 

1090 

19:35:58.73 

-1-78:21:24.3 

0.50 

0 

Y 


* 

0.54 

Y 

AG 

63 

090709A 

253 

19:19:42.64 

-1-60:43:39.3 

0.50 

I 

Y 

Y 

64 

1 s+O-5 
r-o_o.7 

N 

host 

16 

090715B 

57 

16:45:21.63 

-1-44:50:21.0 

0.40 

0 

Y 


* 

3.00 

Y 

AG 

65 

090812 

57 

23:32:48.56 

-10:36:17.2 

0.40 

0 

Y 


66 

2.452 

Y 

AG 

67 

090814A 

13 

15:58:26.39 

-1-25:37:52.4 

0.40 

0 

Y 


* 

0.696 

Y 

AG 

68 

090926B 

71 

03:05:13.93 

-39:00:22.2 

1.40 

X 

N 

Y 

1 

1.24 

Y 

AG 

69 

091018 

14 

02:08:44.63 

-57:32:53.8 

0.50 

0 

Y 


* 

0.971 

Y 

AG 

70 

091029 

24 

04:00:42.62 

-55:57:20.0 

0.50 

0 

Y 


* 

2.752 

Y 

AG 

71 

091109A 

16 

20:37:01.82 

-44:09:29.6 

0.40 

0 

Y 


* 

3.076 

Y 

AG 

72 

091127 

84 

02:26:19.89 

-18:57:08.5 

0.55 

u 

Y 


35 

0.490 

Y 

host 

73 

091208B 

32 

01:57:34.10 

-1-16:53:22.6 

0.25 

0 

Y 


* 

1.0633 

Y 

AG 

74 

100305A 

15 

11:13:28.07 

-1-42:24:14.3 

1.00 

X 

N 

Y 

11 


N 



100615A 

50 

11:48:49.34 

-19:28:52.0 

0.70 

X 

N 

Y 

75 

1.398 

N 

host 

76 

100621A 

210 

21:01:13.08 

-51:06:22.5 

0.40 

I 

Y 

Y 

77 

0.542 

Y 

host 

78 

100728B 

17 

02:56:13.46 

-1-00:16:52.1 

0.52 

u 

Y 


35 

2.106 

Y 

AG 

79 

100802A 

36 

00:09:52.38 

-1-47:45:18.8 

0.50 

0 

Y 


80 

<3.1 

N 

host 

* 

100814A 

90 

01:29:53.59 

-17:59:43.5 

0.40 

0 

Y 


* 

1.44 

Y 

AG 

81 

110205A 

170 

10:58:31.10 

-1-67:31:30.5 

0.30 

0 

Y 


82 

2.22 

Y 

AG 

83 

110709B 

94 

10:58:37.11 

-23:27:16.7 

0.70 

R 

N 

Y 

84 

2.09? 

N 

host 

* 

120119A 

170 

08:00:06.94 

-09:04:53.8 

0.30 

0 

Y 

Y 

85 

1.728 

Y 

AG 

85 

120308A 

12 

14:36:20.05 

-1-79:41:12.2 

0.55 

u 

Y 


35 

<3.7 

N 

host 

* 


^ Swift-BAT prompt-emission fluence (15-150 keV) in units 
of 10“^ erg cm“^. From (ISakamoto et al.ll2011l ) if available, 
otherwise from the Swift GRB table. 

^ Best afterglow position (J2000), relative to the 2MASS 
astrometric system. 

Position uncertainty (arcsec), including an estimate of the 
systematic uncertainty. 

First letter of the wavelength at which the position was 
reported: X-ray, UV, Optical, IR, or Radio/millimeter. 

® Whether or not a variable optical afterglow was reported. 
“IR” indicates a NIR afterglow was reported but not an op¬ 
tical (A < l^m) afterglow. Blank if no observations were 
conducted or reported limits are very shallow. 

^ Whether the afterglow shows evidence of being dust- 
obscured and/or “dark”, based on /lox or the optical-NIR 
color. From a variety of sources including iPerlev et al.l 
(1201.11 1. [Zafar et aTI (|2012l l. [Greiner et al.l (120111 1. 

® Whether or not a redshift estimate was promptly publicly 
available. (If “N”, the redshift was measured by late-time 
observations of the host galaxy or only a limit could be pro¬ 
vided.) 

Source of the redshift measurement/limit: observations of 
the optical afterglow (AG), of the host galaxy (host), or of 
the X-ray afterglow (Ah,x). 

‘ References for the afterglow position (middle column) and 
redshift (right column), given below. 


Sou rces as follows— *: This wo rk. 1: [Evans et al.l 
[20091 2: Malesani, priv. com m., 3: [Berger et al.l [20051 


6: [Fvnbo et al.[[2009bl 

7: Folev et al.l[2005i 8: Malesani 

et al., in prep. 9: ICenkdl2005. 10: 

Jakobsson et al.ll2006L 

II: [Butlei][2007l 12: Hiorth et al.[ 

2012L 13: Kawai et al. 

20061 14: Jakobsson et al. [200f 

. 15: IKriihler et al. 

201 2 I 16: IPerlev et al.ll20L1l 17: 

Guidorzi et al.l 

20061 

18: [Mundell et al.l 120061 19: IMirabal & Haloernl 

20061 

20: iKriihler et al. 120151 21: lPricdl2006. 22: iCenko et al. 

2006al 23: [Levan et al.[ [20071 24: [Ruiz- Velasco et al. 

20071 25: Cenko & Fox [2006[. 26: Perlev et al. 

2009cl 27: IRol et al.l 

20071 28: 

[Jaunsen et al.[ 

2008, 

29: [Chen et al. 2007 

30: iJakobsson et al. 

2007, 

31: IBloom et alJ [^OOtT 

32: ICastro-Tirado et al.l 

2007, 


IBlooml[2007l 34: [pSle 


33 

Table 36: 
[200M 38: 


Svensson et al 


iv et al.[ 

\ [2014 


[mol 35: Swift GRB 
37: [Ghornock et al.[ 


40 

42 

44 

46 


Kriihler, priv. comm. 


39: iHashimoto et al. 

2010, 

41: iRinner & Kngell 

2008, 

43: [Kriihler et al.[ 

2008, 

, 45: |D’Elia et al.l 

1008a, 

47: [D’Avanzo et al. 

l2^. 
















































































































Perley et al. 


22 


48 


Id ’E lia et al.l l2008bl 49: iBerger fc Rauchl 1200 


50: 

Salvaterra et alJ 20121 51: ICucchiara et al. 

|2008b[ 

52: 

Chornock et al.l l2009bl 53: IChornock et al. 

l2009dL 

54: 

Castro-Tirado et al. 20091 55: Berger & Fo 

d 2ooa 

56: 

iPavlenko et al.l l2009l 57: IChornock et alJ 

[2009a[ 

58: IChornock et al.l l2009d 59: Ide Ugarte Postigo et al. 



jTh6nfi^t alJ 


60: 

62: 

64 :^ _ 

66: ICenki 
68 : 

70: IChen et~ 


Goldoni et al 

iMorgan et al.l 


200 ^ 


12009 


67 : 

iJakobsson et alJ 


■201 

200 ^ 


61: 

63: 


iThone et al 


ICenko et al 


65: iWiersema et^ 


de Ugarte Postigo et al 


200^ 72: lR,a,n e.t al.l [20T(1 




|2009l 69; _ 

71: IChornock et al 


Fvnbo et al. 


74 

76 

78 

80 


Perlev et alJ 12009 


Krtihler et al.l 201 


73: ICucchiara et al l [2001 
75: iButler et al. [20 
77: lUpdike et alJ 


Milvang-Jensen et al.l 201 OL 79 : Flores et alJ 


de Ugarte Postigo et alJ l2009dL 81: lO’Meara et aP 



120101 82? iMundell et alJ I2011L 83: fCenko et alJ I2011L 


84: IZanderer fc BergeHl201lL 85: iMorgan et alJl20l¥ 
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The Swift GRB Host Galaxy Legacy Survey 
TABLE 5 

Host Galaxy Photometry Used in Photometric Redshifts and Upper 

Limits 


GRB 

Filter 

mag^ 

unc. 

Instrument 

050128 

R 

24.84 

0.10 

VLT/FORSl 

050803 

g 

>27.74 


GTC/OSIRIS 


R 

26.32 

0.22 

VLT/FORS2 


i 

26.45 

0.50 

GTC/OSIRIS 


F160W 

25.79 

0.18 

HST/WFC3-IR 


3.6 

25.49 

0.42 

Spitzer/IRAC 


4.5 

>25.00 


Spitzer/IRAC 

050922B 

g 

27.63 

0.50 

Keck/LRIS 


R 

26.44 

0.22 

Keck/LRIS 


i 

25.25 

0.14 

GTC/OSIRIS 


z 

25.09 

0.34 

GTC/OSIRIS 


Ks 

>22.17 


VLT/ISAAC 


3.6 

24.59 

0.22 

Spitzer/IRAC 


4.5 

24.80 

0.42 

Spitzer/IRAC 

060204B 

u 

25.33 

0.15 

Keck/LRIS 


B 

24.93 

0.06 

Keck/LRIS 


g 

24.49 

0.03 

Keck/LRIS 


R 

24.12 

0.07 

Keck/LRIS 


i 

23.86 

0.07 

Keck/LRIS 


z 

24.20 

0.18 

Keck/LRIS 


Y 

23.22 

0.16 

Keck/MOSFIRE 


J 

22.46 

0.40 

Keck/MOSFIRE 


Ks 

20.35 

0.40 

Keck/MOSFIRE 


3.6 

22.74 

0.05 

Spitzer/IRAC 


4.5 

22.00 

0.30 

Spitzer/IRAC 

070808 

g 

>27.50 


Keck/LRIS 


R 

26.71 

0.33 

VLT/FORS2 


Ks 

21.77 

0.37 

VLT/ISAAC 


3.6 

23.57 

0.10 

Spitzer/IRAC 


4.5 

23.80 

0.20 

Spitzer/IRAC 

080319A 

u 

25.30 

0.13 

Keck/LRIS 


B 

25.02 

0.10 

Keck/LRIS 


g 

24.50 

0.04 

Keck/LRIS 


R 

24.06 

0.12 

Keck/LRIS 


i 

24.24 

0.08 

Keck/LRIS 


z 

24.26 

0.12 

Keck/LRIS 


Y 

23.61 

0.17 

Keck/MOSFIRE 


J 

22.60 

0.12 

Keck/MOSFIRE 


H 

21.64 

0.12 

Keck/MOSFIRE 


Ks 

21.05 

0.13 

Keck/MOSFIRE 


3.6 

22.63 

0.03 

Spitzer/IRAC 


4.5 

22.60 

0.04 

Spitzer/IRAC 

081128 

B 

26.14 

0.24 

Keck/LRIS 


g 

25.80 

0.15 

Keck/LRIS 


V 

25.37 

0.13 

Keck/LRIS 


i 

25.24 

0.19 

Keck/LRIS 


z 

24.97 

0.33 

Keck/LRIS 


3.6 

25.20 

0.31 

Spitzer/IRAC 

100802A 

B 

26.44 

0.16 

Keck/LRIS 


R 

25.23 

0.25 

Keck/LRIS 


i 

26.23 

0.32 

Keck/LRIS 


3.6 

25.47 

0.43 

Spitzer/IRAC 

120308A 

B 

>26.67 


Keck/LRIS 


g 

26.33 

0.17 

Keck/LRIS 


R 

25.77 

0.22 

Keck/LRIS 


I 

24.46 

0.23 

Keck/LRIS 


3.6 

24.98 

0.27 

Spitzer/IRAC 


In the SDSS (for ugriz), Vega (BRIY), or 2MASS (JHKs) 
magnitude systems; AB magnitudes are used for the IRAC 
(3.6 and 4.5) filters. Measurements are uncorrected for fore¬ 
ground extinction. 
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